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Abstract 

Consider an iV-dimensional Markov chain obtained from TV one- 
dimensional random walks by Doob /t-transform with the g-Vandermonde 
determinant. We prove that as N becomes large, these Markov chains 
converge to an infinite-dimensional Feller Markov process. The dynam- 
ical correlation functions of the limit process are determinantal with an 
explicit correlation kernel. 

The key idea is to identify random point processes on Z with q-Gibbs 
measures on Gelfand-Tsetlin schemes and construct Markov processes on 
the latter space. 

Independently, we analyze the large time behavior of PushASEP with 
finitely many particles and particle-dependent jump rates (it arises as a 
marginal of our dynamics on Gelfand-Tsetlin schemes). The asymptotics 
is given by a product of a marginal of the GUE-minor process and geo- 
metric distributions. 
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1 Introduction 

Let X(t) = (Xi(t), . . . , .Xjv(£)) £ Z N be N > 1 independent rate one Poisson 
processes started at X(0) = (0, 1, . . . , N— 1) and conditioned to finish at X(T) = 
Y while taking mutually distinct values for all < t < T. Thus, {Xi(t)}i>i 
form N nonintersecting paths in Z x [0, T]. 

A result of [25], based on a classical theorem of Karlin and McGregor 
|25j . says that as T — > oo with Y being asymptotically linear, Y ~ T£ with 
a collection of asymptotic speeds £ = (£i < • • • < £jv) € K>o, the process 
(X(t),Q < t < T) has a limit (X(t),t > 0), which is a homogeneous Markov 
process on 7L N with initial condition A(0) and transition probabilities over time 
interval t given by 

P t ((xi, . . .,x N ) -> (yx, . . .,y N )) = const ■ Z N< f \ det 

with Sjv(m) = det [C" 3 ] ^ [^] The process is the Doob /i-transform of 
independent Poisson processes with respect to the harmonic function h = Sjy- 

The case of equal speeds £j = const is special. In that case, the distribution 
of X(t) for a fixed t > is known as the Charlier orthogonal polynomial en- 
semble, which is the basic discrete probabilistic model of random matrix type. 
If one considers its limiting behavior as N —¥ oo. in different parts of the state 
space via suitable scaling limits one uncovers discrete sine, sine, and Airy de- 
tcrminantal point processes which play a fundamental role in Random Matrix 
Theory, cf. PI, HI- 

The procedure of passing to a limit by increasing the number of particles 
and scaling the space appropriately in models of random matrix type is very 
well developed. In many cases such a limit can also be successfully performed 
for joint distributions at finitely many time moments, if the original model 

1 If some of £j's coincide, for the formula to make sense one needs to perform a limit 
transition from the case of different speeds. 



(Vi-XiV. 
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undergoes a Markov dynamics. The most common approach is to control the 
limiting behavior of local correlation functions for the model at hand. 

It is much less clear what happens to the Markovian structure of the dy- 
namics under such limit transitions. While there are no a priori reasons for the 
Markov property to be preserved, it is a common belief that in many cases it 
survives the limit. However, providing a proof for such a statement seems to be 
difficult. 

One goal of the present paper is to prove that if £j 's form a geometric pro- 
gression with ratio q^ 1 > 1, the limit X(t) of X(t) (and other similar processes, 
see below) as N — > oo can be viewed as a Feller Markov process on all point con- 
figurations in Z> started from the densely packed initial condition. Note that 
no scaling is needed in the limiting procedure. As all Feller processes, our limit- 
ing process can be started from any initial condition, and it has a modification 
with cadlag sample paths. 

We also show that the dynamical correlation functions of X(t) started from 
the packed initial condition are detcrminantal with an explicit correlation kernel, 
and they are the limits of the corresponding correlation functions for X(t). 

The process can be interpreted as a restriction of a 'local' Markov 

process on the infinite Gelfand-Tsetlin schemes that falls into a class of such 
processes constructed in [5]. This interpretation implies, in particular, that 
the distribution of the leftmost particle of X(t) coincides with the asymptotic 
displacement of the mth particle, as m — > oo, in the totally asymmetric simple 
exclusion process (TASEP) with jump rates depending on particles, particle j 
has rate and step initial condition. [^] Inspired by this connection, we derive 
large time asymptotics for TASEP, and more general PushASEP of [4_, with 
finitely many particles that have arbitrary jumps rates. 

In the situation when the correlation functions have a determinantal struc- 
ture, the question of Markovianity for systems with infinite many particles has 
been previously addressed in [53], [25], [30] for the sine process, in [27] for 
the Airy process, in [21] for nonintersecting Bessel processes, in (3j5] for the 
Whittaker process, and in [9 for a process arising in harmonic analysis on the 
infinite-dimensional unitary group. 

Our work seems to be the first one that deals with the discrete limiting state 
space. Structurally, we adopt the approach of [5]: We use the fact (proved 
earlier by one of us, see [T7J Theorem 1.1]) that infinite point configurations in 
Z with finitely many particles to the left of the origin can be identified with 
ergodic q-Gibbs measures on infinite Gelfand-Tsetlin schemes. We further show 
that the Markov processes X{t) for different iV's are consistent with respect to 
natural projections from the A-particle space to the (N — l)-particle one; the 
projections are uniquely determined by the q-Gibbs property. Together with 
certain (nontrivial) estimates, this leads to the existence of the limiting Feller 
Markov process. One interesting feature of our construction is that we need to 
add Gelfand-Tsetlin schemes with infinite entries in order to make the space of 
the ergodic g-Gibbs measures locally compact. 

It is worth noting what happens in the limit q — > 1. The space of ergodic 
1-Gibbs measures has countably many continuous parameters (as opposed to 
discrete ones for q < 1), and it is naturally isomorphic to the space of indecom- 
posable characters of the infinite-dimensional unitary group, and to the space 

2 The step initial condition means that we place the particle j at location —j at t = 0. 



of totally positive doubly infinite Toeplitz matrices (see e.g. [17], [33], [35] and 
references therein.) The q — > 1 limit of our Feller Markov process ends up being 
a deterministic (in fact, linear) flow on this space. 

It is plausible that our results on the existence of limit Markov processes can 
be extended to the case of arbitrary {£,j}j>i which grow sufficiently fast, but at 
the moment our proofs of a number of estimates rely on the fact that {£j}j>i 
form a geometric progression. 

We now proceed to a more detailed description of our work. 

1.1 Extended Gelfand-Tsetlin schemes and g-Gibbs mea- 
sures 

Following [48 , for N > 1 we define a signature of length TV as an JV-tuple of 
nonincreasing integers A = (Ai > • • • > Xn)i and we denote by GTV the set 
of all such signatures. For A G GT^r and v G GTjy+i, we say that A -< v if 
v j+i < Xj < v j f° r a U meaningful values of indices. We agree that GTo is a 
singleton 0, and -< v for any v € GUV 

A Gelfand-Tsetlin scheme of order M £ {0, 1, 2, . . . } U {oo} is a length M 
sequence 

A<°> -< A« -< A< 2 > -< A«eGT> 

Equivalently, such a sequence can be viewed as an array of numbers {Ap^} 

satisfying the inequalities A^ 1 < x\ < \+ . An interpretation of Gelfand- 
Tsetlin schemes in terms of lozenge tilings or stepped surfaces can be found in 
[5, Section 1.1]. 

Define an extended signature of length N as an iV-tuple A = (Ai > • • • > Xn) 
with Xj 6 Z = ZU {+oo}R and define an extended Gelfand-Tsetlin scheme as a 
sequence of extended signatures A^ -< A^ -<! A^ 2 ) -(..., A^ € GTj, with the 
condition that the number of infinite coordinates rrij of A") has the property 
that if rrij > then = TOj — 1. In other words, the number of finite 

coordinates grows from level to level until the first infinite coordinate appears, 
on all further levels the number of finite coordinates remains the same. 

Fix q € (0,1). We say that a probability measure on (extended) Gelfand- 
Tsetlin schemes is q-Gibbs if for any N > 2, the conditional distribution of 
^(o) _^ ^(i) _/ . . . _^ / \( JV - 1 ) given A^ w ^ has weights that are proportional to 
gl A< 'l+l A< 'l -1 h l A< 'l, where |A| denotes the sum of all finite coordinates of A. 

There are at least two motivations for our interest in g-Gibbs measures. 
On one hand, if one views Gelfand-Tsetlin schemes as stepped surfaces, the 
conditional distribution above has weights proportional to q volume j where the 
exponent is the volume underneath the surface, and measures of such type have 
various interesting features, cf. [IS], [TT], [35], [B]. On the other hand, the notion 
of 1-Gibbsianness naturally arises in the representation theory of the infinite- 
dimensional unitary group, and we hope that g-Gibbs measures will arise in 
suitably defined representations of inductive limits of quantum groups. 

In |17l Theorem 1.1], one of the authors obtained a description of the convex 
set of q-Gibbs measures on (non-extended infinite) Gelfand-Tsetlin schemes. 
The result says that the space of such g-Gibbs measures is isomorphic to the 

3 We assume that +oo > k for every integer k. 
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space of Borel probability measures on the space 

TV = {vi < v 2 < . ■ ■ , Vi € Z} 

equipped with the topology of coordinate-wise convergence. Moreover, for any 
extreme g-Gibbs measure fi (that corresponds to the delta measure at a point 
v € TV), the value of i/fc, k > 1, is the almost sure limit of A^2 fc+1 as N — > oo, 
with {\^}j>o distributed according to u. 

In this paper we generalize this statement and prove that the space of q- 
Gibbs measures on extended infinite Gelfand-Tsetlin schemes is isomorphic to 
the space of Borel probability measures on the space 

77 = {u 1 <v 2 <..., i/,eZ} 

equipped with the topology of coordinate- wise convergence. Observe that TV is 
a locally compact topological space, and it should be viewed as a completion of 
TV. 

Under the one-to-one correspondence 

W < V2 < v z < . . . } i — > {vi < v 2 + 1 < v 3 + 2 < . . . }, (1.1) 

Af turns into the set of all infinite subsets of Z that are bounded from below, 
and TV can be viewed as the set of all (finite and infinite) subsets of Z bounded 
from below. 



1.2 Consistent N- dimensional random walks 

Let g(x) be a finite product of elementary factors of the form 

(1 - ax- 1 )- 1 , < a < 1, (1 + /3x ±1 ), /3 > 0; exp^a^ 1 ), 7 > 0, 

and g(x) — J2k£i,9k xk be its Laurent series in {x 6 C : 1 < < 00}. This 
means that {<7fe} fc6Z is, up to a constant, a convolution of geometric, Bernoulli, 
and Poisson distributions. 

For any such g(x), define a transition probability on GTjv by 

P N (X ->■ n;g(x)) = const ■ det N \g^-x J+J ] J[ - , (1.2) 

' "' l<i<j<N q q 

and extend it to a transition probability Pat(A — > /i;g(x)) on GT^r by applying 
Pj-( • ; g(x)) with appropriate k < N to all k finite coordinates of extended 
signatures. 

We prove that if one starts with a g-Gibbs measure on extended Gelfand- 
Tsetlin schemes and applies Pn( • , g{%)) to its projections on GT^r, N > 1, then 
the resulting distributions are also projections of a new g-Gibbs measure. This 
gives rise to a Markov transition kernel P ao (g(x)) on TV or, via (1.1 1, on point 
configurations in Z that do not have —00 as an accumulation point. 



1.3 Correlation functions 

Let g (x), gi(x), g2(x), . . . be a sequence of functions as above, and let 2(t), 
t = 0, 1, 2, ... , be the discrete time Markov process on TV with initial condi- 
tion given by the delta measure at = (0,0,...) and transition probabilities 
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Pao(go(x)), Poo(gi(x)), .... One easily sees that the process is in fact supported 
by Af. 

Define the nth correlation function p n of Z(t) at (xi, ti), . . . , (x n , t n ) as the 
probability that Z(tj), viewed via (1.1 1 as a subset of Z, contains Xj for every 
j = 1, . . . , n. We prove that these correlation functions are determinantal: For 
any n > 1, tj € Z>o, Xj € Z, we have 

p n {xi,ti;x 2 ,t 2 ;...;x n ,t n )= det [if (zj, i<; Xj, £,•)], 

— l,...,n 

where 



^(Xi,^;^,^) = -l-j ^™^ \\ g t ( W )l tl 



dw d> dz 



>t 2 

t=t 2 

IltLo 1 9t(w) TT 1 - wq j z X2 



TT- 

(2ttz)2 J c J c , ri^" 1 ft ( z ) -y, 1 - ^ w-z' 

C is a positively oriented contour that includes only the pole of the integrand, 
and C goes from +ioo to —ioo between C and point 1. 



1.4 The Feller property 

A Markov transition kernel on a locally compact state space is said to have the 
Feller property if its natural action on functions on the state space preserves 
the space Co of continuous functions that converge to zero at infinity, cf. [2] . 

We prove that for any function g(x) as above, the kernel P oc (g(x)) on AT 
enjoys the Feller property. 

We also prove that for any 7 + ,7~ > and any probability measure on Af, 
the continuous time Markov process on Af with initial condition fi and transition 
probabilities 

P 00 (exp(i(7 + x + 7 -/x))), t>0, 

is Feller. In addition to the Feller property of the transition kernels, this also 
requires their strong continuity in t when they are viewed as an operator semi- 
group in C {Af). 

Note that such a statement would have made no sense for non-extended 
Gelfand-Tsetlin schemes because Af is not locally compact. This is the place 
where the extended theory becomes a necessity. 



1.5 Asymptotics of PushASEP 

Fix parameters Cij • • • i Civ > 0, a, b > 0, and assume that at least one of the 
numbers a and b does not vanish. Consider N particles in Z located at different 
sites and enumerated from left to right. The particle number n, 1 < n < N, has 
two exponential clocks - the "right clock" of rate a£„ and the "left clock" of rate 
b/Cn- When the right clock of particle number n rings, one checks whether the 
position to the right of the particle is empty. If yes, then the particle jumps to 
the right by 1, otherwise it stays put. When the left clock of particle number n 
rings, it jumps to the left by 1 and pushes the (maybe empty) block of particles 
sitting next to it. All 2N clocks are independent. 
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This interacting particle system was introduced in [4 under the name of 
PushASEP. It interpolates between well-known TASEP and long-range TASEP, 

cf. 021. 

Results of Section 2.7] imply that if Ci = Q 1 ~ N ,C2 = q 2 ~ N , ■ ■ ■ , (n = 
1, and the PushASEP is started from the initial configuration (1 — TV, 2 — 
N, ... ,—1,0), then at time t it describes the behavior of the coordinates 
(A^° +1-N, X^S^ + 2 - TV, ... , A^) of the random infinite Gelfand-Tsetlin 
scheme distributed according to the q-Gibbs measure corresponding to the dis- 
tribution at the same time t of the Markov process on Af with transitional prob- 
abilities {Poo (exp(t(ax + b/x)))}t>o and the delta-measure at as the initial 
condition. 

We prove, for the TV-particle PushASEP with any jump rates and any de- 
terministic initial condition, and independently of the rest of the paper, that 
at large times the PushASEP particles demonstrate the following asymptotic 
behavior: In each asymptotic cluster, particles with the lowest values of 
fluctuate on \fi scale, and the fluctuations are given by the distribution of the 
smallest eigenvalues in an appropriate GUE-minor process, while faster parti- 
cles remain at O(l) distances from the blocking slower ones, with geometrically 
distributed distances between neighbors. The GUE-governed fluctuations and 
the geometric distributions are asymptotically independent. 

Note that for TASEP, which arises when 6 = 0, Ci = C2 = • • ■ = (n, the 
relation of large time asymptotics with marginals of the GUE-minor process is 
well known, see [2J. 



1.6 Organization of the paper 

In Section [2] we discuss TV-dimensional random walks with transition proba- 
bilities P/v(A — ► fi;g(x)). In Section [3] we prove that these random walks are 
consistent for different values of TV. In Section [3] we explain how the family of 
consistent transition probabilities define the transition kernel -Poo (A —> fx;g(x)) 
on Af and describe basic properties of this kernel. The properties of the corre- 
sponding Markov processes are studied in Sections [5] and [6j 

The main results of Section [5] are Theorem |5.1| where we prove that finite- 
dimensional distributions of the constructed processes are TV — >• 00 limits of 
those for the TV-dimensinal processes; and Theorem |5.6| where we prove the 
corresponding result for the correlation functions. In Section |5.3| we describe 
the connection between processes with transition probabilities Poo (A — ► fi; g{x)) 
and stochastic dynamics of [5]. 

Sectio n [6 is devoted to the study of extended Gelfand-Tsetlin schemes. In 
Theorem |6.3 we prove that the space of g-Gibbs measures on extended infinite 
Gelfand-Tsetlin schemes is is omor phic to th e space of Borel probability measures 
on the space Af. Theorems 6.15 and 6.17 contain the proofs of the results on 



Feller properties mentioned above. This is the most technical part of the paper. 

Finally, in Section [7] we prove Theorem 7^2 describing the asymptotic behav- 
ior of PushASEP with finitely many particles. 
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2 N— dimensional random walks 



In this section we introduce a family of Markov chains on the set of iV-point 
configurations in Z. In a variety of situations they can viewed as collections of N 
independent identically distributed random walks conditioned not to intersect 
and to have prescribed asymptotic speeds at large times. 

The state space of our processes is the set M^v of iV-point configurations in 



= {<z(l) < a(2) < • • • < a(N) \ a(i) e Z}. 



We identify elements of Mjv with weakly decreasing sequences of integers, which 
we call signatures. The set of all signatures of size N is denoted by GT/v, 

<GT N = {Ai > ••• > X N | Xi eZ}. 

We use the following correspondence between the elements of Mjv and GT^: 

a(l) < •• • < a(N) i — > a(N)-N + l > a(N-l)-N + 2 > ■■■> o(2)-l > a(l). 

We also agree that Mo (GTq) consists of a singleton — the empty configuration 
(the empty signature 0). 

'Signature' is the standard term for the label of an irreducible representation 
of the unitary group U(N) over C, and the letters GT stand for 'Gelfand-Tsetlin' 
as in Gelfand-Tsetlin bases of the same representations. Although the material 
of this paper is not directly related to representation theory, we prefer not to 
change the notation of related previous works, cf. |36l \17\. 

In studying probability measures and Markov processes on GT^v we ex- 
tensively use rational Schur functions. These are Laurent polynomials s\ £ 



] indexed by A G GTjy and defined by 



Sx(x!,. . . ,X N ) = - 



rii^^i X j) 

Let us introduce certain generating functions of probability measures on 
GT n that may be viewed as analogues of characteristic functions but that are 
more suitable for our purposes. 

Fix a non-decreasing sequence of positive reals {(i}i=i 2,...' Let Tjv be the 
TV-dimensional torus 

T N = {( Xl ,...,x N ) eC N I W = &}. 

Denote by Tn a class of functions on T/v which can be decomposed as 

f( Xl ,..., XN )= £ c x (f) S ^-- XN) (2.1) 

with ca(/) > and J2\ c ^(f) = !• Note that the latter condition is equivalent 
to /(&,..., frr) = l. 



Lemma 2.1. The series (2.1) converges uniformly on Tjs and its sum 
f{x\, . . . , a;jy) is a real analytic function on T/v- 
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Proof. This fact immediately follows from a simple observation 

s\(xi, . . .,X N ) 



sup 

(a:i,...,xjv)eTw 



sa(£i, ■ • ■ ,&v) 



= 1, 



which, in turn, follows from the combinatorial formula for Schur functions (see 
e.g. Section 1.5 of (32]). □ 

The authors know no intristic definition of the set J-jy and it would be 
interesting to find one. 

We note without proof that Fn is a closed subset of the Banach space of 
continuous function on Tjv, see |17l Proposition 4.11] for a similar fact and |17l 
Section 6.2] for its proof which translates to our case almost literally. 

Let P be a probability measure on GIV- Its Schur generating function is a 
function S(x±, . . . , x^; P) £ J-jv with coefficients c\(S) defined through 

c x (S) = P(X), 

where P(X) stays for the measure of the singleton {A}. Let Cn be the map 
sending probability measures on <GT/v to the corresponding functions in J-jy. 
Clearly this is an isomorphism of convex sets. We agree that J-q contains a 
single function (constant 1) which corresponds to a unique probability measure 
on the singleton GTq. 

Our next goal is to construct a family of stochastic matrices with rows and 
columns enumerated by elements of <GTV. Let Q be one such stochastic matrix. 
Then Q acts on probability measures on G1V 

P^QP. 

We will always identify stochastic matrices with the corresponding operators 
and use the same notations for them. 

Let Q be a bounded linear operator in the Banach space of continuous func- 
tions on Tn such that Q(-^A') C J~n- Then CJ^QCn is a Markovian linear 
operator or a stochastic matrix. 

For a function g(x) on {jf =1 {x E C : \x\ — define an operator Q 9 N via 

N 



f(Xi, . . .,X N ) I ^ f(xi, . . . ,x N ) Y[ ' 



(here we agree that Qq is the identity operator). Clearly, if the function g(x) 



is continuous then Q 9 N is a bounded linear operator in the space of continuous 



functions on T/v- 

Lemma 2.2. Suppose that g(x) can be decomposed into a converging power 
series in annulus K = {x £ C : r < \x\ < R}: 

g{x) = ^c k x k , i£K. 

fcez 

N 



Then for (x%, xn) & K we have 



8\(xi, . . . 1 x N )g(x 1 ) ■ ■ -g{x N ) = V" det 

— ' i,j — l,...,N 



-i—Xj+j 



S fl (xi, . . . ,X N ). 

(2.2) 
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Proof. Straightforward computation. One multiplies both sides of (2.2 1 by 
6.5] 

Denote by g(N, £) the set of functions consisting of 



nj<j( a; j — x j) anc ^ compares the coefficients of the monomials, cf. [381 Lemma 

□ 



x; 

(1 + px), 
(1 - ax)- 1 



-l. 



P > 
7 > 
< a < Q 1 



(1+Px- 1 ), 



/3>0; 
7>0; 

(1-ax -1 ) -1 , < a < £i. 



Also denote fl(oo,£) = f) N q(N,£). In what follows we call g(x) € g(oo,£) an 
elementary admissible function. Let <7(./V, £) denote the set of all finite products 
of the functions of q(N, £) and denote <?(oo,£) = C\n^(^'0- We call g(x) S 
£/(oo,£) an admissible function. 

Proposition 2.3. //ye G(N,£), then 

Q%{Fn) C .Fat. 



Proof. For g(a;) € G(N,£) we can use Lemma 2.2 Furthermore, it is known 
that in these cases all determinants in the decomposition (2.2 1 are non-negative 

(see m, [SI, HO], @7], (36]). 

Then we proceed as follows: Take a function / € .Fjv, decompose it into a 
sum of Schur polynomials. Then by Lemma |2.2| we obtain a double sum for 
Q 9 N {f). Changing the order of summation we see that Q%{f) £ J~n- D 



Set 



P N (g(x))=£^oQ 9 N ojC 



-A' 



and let Pat(A — > fi;g(x)) be the matrix element of the corresponding stochastic 
matrix. 



Proposition 2.4. Let 



g{x) = ^2 c k xk 



be a decomposition of g{x) into a power series converging for all x such that 
mini £i < \x\ < maxi £j, i = 1, 2 . . . , N. Then 



P N (X^ ti;g{x)) = [fl^Tj 



det 

ij=l,...,N 



Proof. This is an immediate consequence of Lemma |2. 2 1 and definitions. 



(2.3) 



□ 



Remark 1. For g(x) = (1 — ax) 1 the determinants in the proposition 
above can be explicitly evaluated: 



det 

i,j = l,...,N 



| a E£ 1 (w-A«) j if < A . < M . ; x < j < Nj 

[ 0, otherwise. 

(2.4) 

(The condition /iq < Ai above is empty.) 
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For g(x) = 1 + (3x, the evaluation takes the form 

f/3Ef =1 (w-AO, if ^ _ Xi e {0, 1}, 1 < i < N, 
] 0, otherwise. 



(let 

i,j=l,...,N 



Similar formulas exist for g{x) = (1 — a 1 x) 1 and g(x) = 1 + fix 1 with A 
and /i interchanged. 

Remark 2. When £ is a geometric progression, the Schur function 
sa(£i, ■ ■ ■ ,£n) can be evaluated as follows (see e.g. [32, Example 3.1]) 

S X (1, . . . , q 1 -") = q -((N-l)X 1 + (N-2)X 2+ - + X N ^) -Q ' ^ ' '] ' ' . (2 .5) 



Then the formula for the transition probability (2.3 1 turns into (1.2 1. 



Remark 3. If & = q 1 ' 1 and g(x) = (1 + fix* 1 ), then one can formally send 



TV — > oo in formulas (1.2), ( 2 .3 1 and obtain well-defined transition probabilities; 
while for g(x) = exp(7a;) such formal limit transition does not lead to anything 
meaningful. 

Denote by XN, g (t) a discrete time homogeneous Markov chain on CTjv with 
transition probabilities Pat(A — » fi;g(x)), started from the delta-measure at the 
zero signature = (0 > > • • • > 0). 

Also let 3V l7+l7 _ (t) be the continuous time homogeneous Markov chain on 
GT^v with transition probabilities Pjy (A — > /x; exp(f(7 + x + 7_:r~ 1 ))) , started 
from delta-measure on zero signature 0. (Clearly the corresponding stochastic 
matrices form a semigroup.) 

A number of such Markov chains with various g(x), 7+, 7_, have independent 
probabilistic interpretations. Let us list some of them. 

• For N = 1 and £i = 1, Xi t i+p x (t) and X N i + p x -i (t) are simple Bernoulli 
random walks with jump probability /3(1 + /3) _1 and particle jumping to 
the right and to the left, respectively; X lx ±i(t) is the deterministic shift 
of the particle to the right (left) by 1; X!in_ ax \-i(t) and X N (i_ ax -i^-i (i) 
are random walks with geometrical jumps; 3^i.7 + ,o(^) is the Poisson process 
of intensity 7 + and J^i,7 + ,7_ (t) is the homogeneous birth and death process 
on 1. 

• For any N > 1 and an arbitrary sequence £, it is proved in [29J that 
3^v,7 + ,oW can be viewed as N independent rate 1 Poisson processes con- 
ditioned never to collide and to have asymptotic speeds of particles given 
by Similar interpretations could be given for X^ g (t). 

• If £i = £2 = • • • = £w = 1, then A'jv,i+/3x(i) can be obtained as a limit of 
uniformly distributed 3d Young diagrams in a x b x c box (see [H5], |20| . 
[22_, [16 ) with a = N and b, c — > 00 in such a way that c/b — > /3. 

• The connection to exclusion processes is explained in Section |5 . 3| below . 



Proposition 2.4 implies that the one-dimensional distribution of X^ tg (t) at 



a given time to is a (possibly, two-sided) Schur measure, cf . [3] , 



Prob(AW ig (i ) = A) 




to 



-t+J 



sa(£i,---,£ 
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where are the coefficients of the Laurent series for (<?(#)) °. 

If we view X^ g (t ) as a point configuration in Z, i.e. as an element of M^r, 
then we may speak about its correlation functions: 

p n (xi, ...,x n ) = Prob(a;i G <*V, ff (io), • ■ ■ ,i n € X N<g (t )). 

As shown in [7] , the correlation functions of the two-sided Schur measures 
have a determinantal form (for the one-sided Schur measure this was proved 
earlier, see pH]) 

p n (xx,...,x n ) = det [K(xi,Xj)], 

l,j=l,...n 

with correlation kernel K given by a double contour integral 



K(x 1 ;x 2 ) = -^fd W f dz ^ J] f x 1+ i — ■ ( 2 - 6 ) 

(2my J c J c < \g{z) J fj- (1 - z/&) w*^ 1 w - z 



where the (positively oriented) contour C includes only the pole at and C 
includes only the poles at £j. A similar formula exists for dynamical corre- 



lation functions (see P3 Section 2.8] and Section 5.2 below), describing finite- 
dimensional distributions of XN. g {t) and also for D^i,7 + ,7_ (£)■ 

Note that if £j's grow fast enough as i —¥ oo, then one may formally pass to 



the limit N oo in (2.6). Thus it is natural to expect that there is an infinite 
dimensional Markov process which is a N — > oo limit of processes Xff <g (t) (or 
3^jV,7 + ,7_ (*))■ ^ ne S oa ^ °^ * n * s P a P er i s to define this limiting process rigorously 
and to show that its finite-dimensional distributions are indeed given by limits 
of KE). 



3 Commutation relations 

In this section we show that transition probabilities P/v(A — > /j,;g(x)) are in a 
certain way consistent for various N. 

We start by introducing stochastic matrices with rows and columns indexed 
by elements of GTjv and GTjy-i, N > 1. In other words, we want to define 
transition probabilities from GTTV to GTjv-i- As above, it is convenient to pass 
from stochastic matrices to maps between spaces of Schur generating functions 
of probability measures. Thus, we want to introduce a map 

Proposition 3.1. The specialization map 

P N '■ f( x i>- ■ -,xn) -> f(xi, . ■ .,x N -i,£ N ) 

is a bounded linear operator between appropriate spaces of continuous functions, 
and P^^n) C 7jv_i. 

Proof. The fact that this is a bounded linear operator is straightforward. Using 
well-known branching rules for Schur functions (see e.g. |32| ) we see that: 

pi f s\(xi,...,x N ) \ _ s^jxi, . . . ,xat_i) |a|-|mI s P (Ci, ■ ■ -,£n-i) 
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where fi -< A means the following interlacing condition: 

Ai > Hi > A 2 > • • • > fiN-i > Ajv, (3.1) 
and |A| = Ai + • •• + Ajy, |a*I = Ml + "" + MiV-i- Since all the coefficients 



\\\-\h\ s,i(€u . . . ,£n-i) 



are positive, we immediately conclude that P^J-jv) C J-jv-i- O 

Let us denote by Pjy a stochastic matrix of transition probabilities corre- 
sponding to P^r, i.e. Pj^ = C^_i ° Pjy o Ln. We call this matrix a stochastic 
link between levels N and N — 1 . 

Using the definition of Schur functions we conclude that the matrix elements 

n( 



P^(A — > ji) are given by the following formula: 



0, otherwise. 

(3.2) 

Note that if £j's form a geometric progression then the determinants in (3.2 1 



turn into g-Vandermondc determinants, cf. (2.5 



Proposition 3.2. Matrices of transition probabilities Pjv(- ; g(%)) commute 
with links Pjy. 

Proof. This is equivalent to commutativity relations between maps Q 9 N and Ph, 
which is straightforward. □ 

For our further constructions it is necessary to extend the space GTjv and 
the definition of links Pj^ . 

The extended level N, GT^ consists of all sequences Ai ^ A 2 2^ ■ * * ^ A^y, 
where Aj 6 ZU {+oo}. We identify every such sequence with a shorter sequence 
Mi > /i2 > • • • > /^fc with /Zj = Ai+Ar_fc and k being the smallest integer such 
that Ajv-fc = +oo, and with the corresponding fc-point configuration in Z. 

We have 

GlF/v = GT^ U GT% ] U GT^ U • • • U GT^' , (3.3) 



where 



GT^ — {A : Ai — • • • — \N-k — +oo, \ N -k+i, . . . , Ajy G Z}, 



in particular, GT^ consists of a single signature with all infinite coordinates. 
It is convenient to use the obvious identification 

GT^ } ~ GT^ ~ • • • ~ GT<: fe) = GT fc . 

In order to define the extended matrix of transition probabilities Pj\r(A — > 
A € GUV, fi e GT^r-i we first introduce for any k < N an auxiliary stochastic 
matrices Q 1 ^ with rows and columns indexed by elements of GT^ by = Q? 
with g — 1/(1 — ^^x). And Q° N is the unique lxl stochastic matrix. 
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Now we are ready to d efin e P j^(A — > fx). This matrix has a block structure 
with respect to splittings (3.3 1 on levels N and N — 1. For A e GT^' 

(A H. M ) = |^ (A ^ M) ' if ^ G GT N-l\ 
[ 0, otherwise. 

For A E GTP with k < N, we define 

5^ , A*), if AiGGT^j, 



0, otherwise. 



Let us also extend the definition of stochastic matrices Pjy( ■ ; g(x)) to larger 
matrices Pn( - ; g(x)) with rows and columns indexed by elements of GT/v- 
These matrices also have a block structure with respect to (3.3). We define 



P N (X -> fi;g(x)) 



' P k (-;g(x)), if A E GT^ and fi € GT^ } , 



0, otherwise, 

where the arguments of P k on the right are suitable truncations of A and ji. 



Proposition 3.3. Matrices of transition probabilities P^{-;g{x)) commute 
with links P^ N . 

Proof. This is immediate from the definitions. □ 



4 Infinite— dimensional dynamics. 

In this section we introduce infinite-dimensional dynamics — the main object 
of study of the present paper. We start from general theorems and then we 
specialize them to the probability measures and Markov chains on the sets 
GT N . 

4.1 General constructions 

For any topological space W we denote by A4(W) the Banach space of signed 
measures on W with total variation norm, and by A4 p (W) the closed convex 
subset of probability measures. 

Suppose that we have a sequence of countable sets ro,ri,. .. , and for any 
N > we have a stochastic matrix A^ +1 with rows enumerated by elements of 
rjv +1 and columns enumerated by elements of Tjsf. We call these matrices links. 
Aj^ +1 induces a linear operator mapping A4(Tn + i) to A4(Tn) and we keep the 
same notation for this operator: 

(A% +1 M)(y)= ]T M(x)A% +1 (x,y), M € M(T N ). 

The projective limit lim A4(IV) with respect to the maps Aj^ +1 is a Banach 
space with norm 

||(M ,Mi,...)|| = sup||Mjv||. 

N 
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Suppose that the sets Tjv are equipped with some topology (which can be, 
in principle, the discrete topology). To avoid pathologies we assume that these 
topologies are metrizable. Then the set ^m.^4(T^) has another natural topology 
which we call the weak topology. The weak topology on limA^(rjv) is the 
minimal topology such that for every N > and every bounded continuous 
function f(w) on Tjv the map 

{M N }^ J2 f(v>)M N (w) 

is continuous. 

We equip )miM.(Tpj) with Borel er-algebra spanned by open sets in the 
norm-topology. One proves that this is the same algebra as Borel c-algebra 
spanned by open sets in the weak topology. (This is the only place where we 
need metrizability. It may be unnecessary here as well but we do not have a 
proof of that.) 

A projective limit hjn A4 p (Tn) is a closed (in both topologies) convex subset 
of ljm A^(rjy). Elements of lim M p (Tn) are called coherent systems. Note that 
if M is a coherent system, then 

||Mo|| = ||Mi||=--. = ||M|| = l. 

Definition 4.1. A topological space Q is a boundary of a sequence {Tn, Aj^ +1 } 
if 

1. There exists a bijective map 

E: M(Q) -> lim M(T N ); 

2. E and E^ 1 are bounded linear operators in the corresponding norms; 

3. E maps A4 P (Q) bijectively onto l^mA4p(T^); 

4- x — > E(S X ) is a bijection between Q and extreme points of the convex set 
of coherent systems, and this bijection is a homeomorphism on its image, 
where we use the restriction of the the weak topology of lim A4(Tn)- 

Remark 1. If Q\ and Q2 are two boundaries, then they are homeomorphic. 

Remark 2. As follows from 4. of the above definition, the boundary can be 
always identified with the set of all extreme coherent systems. However, finding 
a more explicit description of the boundary can be complicated. 

Remark 3. Some authors define the boundary to be the set of all extreme 
coherent systems (see e.g. [5S], [T^], |35]), then 1.-4. become the properties of 
the boundary. Also note that a slightly different definition of the boundary was 
used in jjj]. 

Theorem 4.2. For any sequence {rjv,Aj^ +1 } there exists a boundary Q. 
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Proof. Statements similar to Theorem 4.2 were proved in [T2] and [35]. We use 
|38| as our main reference. 

Let Q be the set of extreme points of the convex set hm A4 p (Tn). We equip 
Q with weak topology inherited from UmA^(rjv). Note that condition 4. of 
Definition |4.1| is satisfied automatically. 

Theorem 9.2 in [35J says that Q is a Borel subset of HmA^ p (rAr). and for 
any M £ lim M p (Tn) there exist a unique probability measure ttm on Q such 
that 

M= [ M«Tr M (dq), 

Jq 

meaning that for any N > and any subset A of Tn we have 

M N (A)= [ M^(A)7T M (dq). 
Jq 



For any tt G set 



E(n) = [ Af q 7r(dq) 

Jq 



If 7r is a signed measure, i.e. tt £ Ai(Q), then there exist %i,TT2 € A^ P (Q) and 
two non-negative numbers c\, c 2 such that 

TT = Cl7Tl — C 2 7r 2 

and ci + c 2 = ||7r||. We define 

E(ir) =c 1 E(n 1 ) - c 2 E{tt 2 ). 

Clearly, E is a linear operator from Ai(Q) to lim A4(Tn) of norm 1. Using 
Theorem 9.2 from [35] we conclude that condition 37of Definition fLll is satisfied 
and, moreover, E is an injection. 

In order to prove that E is a surjection it is enough to show that for any 
M £ luaM{T N ) there exist K,L <E \\mM p {T N ) and k,l > such that M = 
kK — ZL. Without loss of generality we may assume that neither M, nor — M 
are positive measures. Let Mm = — L 1 ^ be a decomposition of Mjv into 
a difference of two positive measures such that ||M;v|| = ||-K"jvH + l^jvll- Set 
K^_ t = Aft^Kft, K^_ 2 = A^ZlKft^ and so on, and similarly for L%. Note 
that for any N > 0, K N monotonically increases as i — > oo. Also < ||M||. 

Hence, there exists a limit 

K% = lim K N . 
In the same way there is a limit 

L% = lim L N . 

i— > oo 

Note that for any i, K N and U N are positive measures and Mrsr = K l N — L N . 
Therefore, similar statements hold for M, and L 1 ^ . Setting k — \\K°°\\, 
I = \\L°°\\ (neither k nor I can vanish because we assumed M is not positive), 
K = K°° /k, L — L°° /I we get the required decomposition of M. 

We have proved that E is a bounded linear operator in Banach spaces which 
is a bijection. Then it follows from Banach Bounded Inverse Theorem that E~ x 
is also a bounded linear operator. □ 
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Next, we want to introduce a way to define a Markov chain on the boundary 
Q of {IV, A^ +1 } using Markov chains on the sets -{Tat}. 

Let Wi, W2 be topological spaces. A function P(u, A) of a point u £ Wi 
and Borel subset A of W2 is a Markov transition kernel if 

1. For any u £ Wi, P(u, • ) is a probability measure on W2, 

2. For any Borel subset A C W2, P( ■ , A) is a measurable function on W\. 

If Wi = W2 = VV, then we say that P(u, A) is a Markov transition kernel on 
W. 

Suppose that for every N > we have a Markov transition kernel on r^r 
given by a stochastic matrix Vn(x y), x,y € Tjv- Assume that these matrices 
commute with links A^ +1 , i.e. for N = 0, 1, 2, . . . we have 

V N A% +1 =A% +1 V N+1 . (4.1) 

Now we define a kernel Voo on Q in the following way. Take a point u £ Q; 



by Theorem 4.2 Dirac 5-measure 5 U corresponds to a certain coherent system 



{Af^r}. Then relations (4.1 1 yield that {Pat Mat} is also a coherent system and, 
consequently, it corresponds to a probability measure on Q. We define • ) 

to be equal to this measure. 

Proposition 4.3. is a Markov transition kernel on Q. 

Proof. Voo(u, • ) is a probability measure on Q by the very definition. Thus, it 
remains to check that Poo( ■ , A) is a measurable function. Theorem |4.2| yields 
that the map u 1— > M^(x) (where = E(S u )n) is continuous and, thus, 
measurable for any N > and any x £ Tjv. Then the map u (VjsrM^)(x) is 
also measurable. Since, the last property holds for any x, we conclude that u m- 
{PnMn} is a measurable map from Q to coherent systems. The definition of 
the boundary implies that the correspondence between probability measures on 
Q and coherent systems is bi-continuous and, therefore, it is bi-measurable. We 
conclude that the correspondence u 1— > T' 00 {u 1 •) is a measurable map between 
Q and M P (Q) (the Borel structure in the latter space corresponds to the total 
variation norm). It remains to note that for any Borel A C Q the map 

Ev A :M(Q)^R, Ev A (P) = P(A) 

is continuous (in the total variation norm topology), thus measurable. There- 
fore, for any Borel A, the map u Voo{u, A) is measurable. □ 

We see that transition kernels Vn naturally define Markov chains on Q. 



However, Theorem 4.2 tells us nothing about the actual state space of these 
Markov chains, and their properties can be very different. Next we state suffi- 
cient conditions for these Markov chains on Q to enjoy the Feller property that 
we now recall. 

Let W be a locally compact topological space with a countable base. Let 
B(W) be the Banach space of real valued bounded measurable functions on W, 
and let Co(W) be the closed subspace of all continuous function tending to zero 
at infinity. In other words, these are continuous function f(w) such that for any 
e > there exist a compact set V C W such that < £ for all w £ W \ V. 

Note that Cq(W) is separable and its Banach dual is A4(W). 
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Let P(u,A), u £ W\, A C W2 be a Markov transition kernel. It induces a 
linear contracting operator P* : B(y\>2) — > B(VVi) via 



(P*f)(x)= / /(«;)P(x,di£;). 

We say that P(it, A) is a FeZfer fcerne/ if P* maps Co(W 2 ) to C (Wi). 

Now let X(t) be a homogeneous continuous time Markov process on W with 
transition probabilities given by a semigroup of kernels Pt{u, A). X{t) is a Feller 
process if the following conditions are satisfied: 

1. P t (u,A) is a Feller kernel, i.e. P ( * preserves Co(W). 

2. For every / £ Co(W) the map / — > P* f is continuous at f = 0. 

Feller processes have certain good properties. For instance, they have a 
modification with cadlag sample paths, they are strongly Markovian, they have 
an infinitesimal generator, see e.g. |14l Section 4.2]. 

Proposition 4.4. Let, as above, Q be the boundary of {T ^ , A^ +1 } , Vn be 
stochastic matrices on Tn commuting with Aj^ +1 , and let Voo be the corre- 
sponding Markov transition kernel on Q. Suppose that 

1. Tjv and Q are locally compact topological spaces with countable bases; 

2. For every N > the function A™(u, A) := E(6 U ) N (A), u £ Q, A CT N is 
a Feller kernel; 

3. Matrices Vn define Feller kernels onT^. 
Then Voo is a Feller kernel on Q. 

Proof. Let h £ C (Q). We need to check that T^(h) £ C (Q). Since C Q (Q) 
is closed and V^ xi (g) is a contraction, it is enough to check this property on a 
dense set of functions h. Let us find a suitable dense set of functions. 

We claim that the union of the sets (k^)* (Cq(T n)) over N > is dense in 
Co(Q). Indeed, A4(Q) is a Banach dual to Co(Q), thus, it is enough to check 
that if 7T £ M{Q) is such that J fdn = for all N and all / £ (A%)*(C (T N )), 
then 7r = 0. The latter property is equivalent (by Fubini's theorem) to the 
following one: For any N and any / £ Co{Tn)- J fdA^(Tr) — 0. But then 
A^(7r) = 0, therefore, E(ir) = and it = 0. The claim is proved. 

Now let h = (A^? )*(/). Then by the definitions 

VUh) = PSo((A??)*(/)) = A%(V* N (f))- 



But T* N (f) £ C (T N ) by the condition 3. of Proposition 4.4 Thus, V^(h) £ 



C (Q) by the condition 2. □ 

Proposition 4.5. LetX^(t) be a homogeneous continuous time Markov process 
onTpj with transition probabilities given by a semigroup of stochastic matrices 
Vt,N- Suppose as above, that Q is the boundary of {Tjv, A^ +1 }, the matrices 
Vt,N commute with links A^ +1 and let Vt,oa be the corresponding semigroup of 
Markov transition kernels on Q. 

If Xjy(t) is a Feller process for every N > and, furthermore, assumptions 
1. and 2. of Proposition \J^4\ are also satisfied, then a Markov process on Q with 
semigroup of transition probabilities Pt jCO o,nd arbitrary initial distribution is 
Feller. 
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Proof. The first defining property of the Feller process is contained in Proposi- 
tion 4.4 As for the second one, since Vt,oo is a contraction, we may check this 



property on a dense subset. If h = (A^)*(/) and / £ Co(Tn), then 

Vt 00 (h) = (A%)*(Vl N (f)). 

And this is a continuous map as a composition of a continuous map and a 
contraction. □ 

4.2 Specialization 

Now let us specialize the general situation by setting = GTjv (with dis- 
crete topology) and Aj^ +1 = P^ +1 , where these matrices implicitly depend 
on the sequence {£tv}- Denote by <2(£i,£2, ■ • • ) the boundary of the sequence 
{GT Nl P^ +1 }. 

Given a n ad missible function g(x) € C/(oo,£), we set Vn '■— Pn(u, A; g(x)). 
Proposition 3.2 yields that these stochastic matrices commute with links Pj^ +1 , 



thus, the constructions of Section |4~T| give as a Markov transition kernel on 
Q(£i:£2, ■ • • ) which we denote Poo(u, A; g(x)). 

As far as the authors know, an explicit description of the set <2(£i,£2, • ■ ■ ) 
is currently known in two special cases. Namely, 

1. £i = £2 = • •• = 1- This case was studied by Voiculescu |44| . Boyer |10j . 
Vershik-Kerov |Uj, Okounkov-Olshanski [36 . It is related to represen- 
tation theory of the infinite-dimensional unitary group and to the classi- 
fication of totally-positive Toeplitz matrices. However, it turns out that 
in this case the Markov operators P oc (g(x)) correspond to deterministic 
dynamics on fi, cf. [3]. 

2- 6v = q^ N , < q < 1. This case was studied in [T7]. As we show 
below, P oc (g(x)) leads to a non-trivial stochastic dynamics. (Note that 
the case q > 1 is essentially the same asO<g<luptoa certain simple 
transformation of spaces G1V.) 

From now on we restrict ourselves to the case £w = q X ~ N ■ N > 1. The 
following results were proven in |17l Theorem 1.1, Proposition 1.2]. 

Theorem 4.6. For £ N = q x ~ N , < q < I, the boundary Q(£i, £2, ■ ■ • ) « s 
homeomorphic to the set M of infinite increasing sequences of integers 

N = {v x < v 2 < Ui GZ} 

with the topology of coordinate-wise convergence. 

Denote by E q the bijective map from (signed) measures on J\f to 
limy\4(GT7v) ; and let £ u — E q (8 u ). Then the coherent system E^, N — 
571)2,..., has the following property: If we view Xn, Xn-i, ■ ■ ■ as random 
variables on the probability space (GTjv,£/v), then for every k > 1, Ajy-fe+i 
converges (in probability) to v% as N — > 00. 

Similarly to G1V, we identify elements of Af with point configurations in Z: 
v\ < U2 < • • • < — > {vi + i — 1, i = 1, 2, . . . }. 
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Note that in this way we get all semiinfinite point configurations in Z; we denote 
the set of such configurations by M^. 

The space Af is not locally compact. This introduces certain technical diffi- 
culties in studying continuous time Markov chains on this space. To avoid these 
difficulties we seek a natural local compactification of Af. 

Let 



Af = Q GT N u Af. 



N=0 

We identify elements of GT n with infinite sequences V\ < . . . such that 
Vi, ... , vn are integers and vn+i = = ■ ■ ■ = +00. Thus, 

A7 = {fi < V2 < ■ ■ ■ , Vi £ Z U {+00}}. 

In the same way we set 

= Moo U M U Mi U . . . 

Clearly, M.^ is a set of all point configurations in Z which have finitely many 
points to the left from zero. There is a natural bijection between Af and M^. 

We equip Af with the following topology. The base consists of the neighbor- 
hoods 

A h k = {v € Af : vx = r/i, . . . , v k = rj k }, rji £ Z, 

and 

Bn,k,i = {v eJV" : V\ = Vu- ■ -i v h = Vk^k+i > £}, % € Z, I £ Z. 

The following proposition is straightforward. 

Proposition 4.7. Topological space Af is locally compact, the natural inclusion 
Af Af is continuous, and its image is dense in Af. 

Now we are ready to define a kernel Poo{ u , A; g(x)) on Af: If u £ GTyv C Af, 
then P oc (u, A; g{x)) is a discrete probability measure concentrated on GTjv with 
weight of a signature A £ GTjv given by Pn{u — > A; 5(2;)); if u £ Af C A/", then 
measure Poo{u, A; g(x)) is concentrated on Af and coincides with Poo(u, A; 5(2;)) 
on it. 

Proposition 4.8. P 0o (u, A; g(x)) is a Markov transition kernel on Af. 

Proof. For every u, Pqo(u, ■ ; g(x)) is a probability measure by the defini- 
tion. The measurability of P ca (- ,A;g(x)) follows from the mcasurability of 
P oc (-,A;g(x)) tmd P N (- -> A;g(x)). □ 

In Section 6 we prove that Poo{u, ■ ; g{x)) is a Feller kern el by identifying Af 
1 

with a boundary of {GTjy,P N } and then using Proposition 



4.4 



5 Description of the limiting processes 

In this section we study finite-dimensional distributions of Markov processes 
that correspond to the Markov kernels Poo(u, A; g(x)). 
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5.1 A general convergence theorem 

Let gk(x), k = 0, 1, 2, . . . , be a sequence of admissible functions, and let Zq be 
an arbitrary probability distribution on Af. 

Denote by Zjy(t), t = 0, 1, 2, ... , a discrete time Markov chain on GT/v with 

initial distribution Zpj(0) = E q (Z$)N and transition probabilities 

Pmb{Z N (t + 1) G A | Z N (t)} = P N (Z N {t), A-g t {x)). 

Also let -Z(i), t = 0, 1,2, . . . , be a discrete time Markov process on Af with 
initial distribution 2(0) = Zq and transition measures 

Prob{Z(t + 1) e A | Z(t)} = Poo(Z(t),A;g t (x)). 

Note that the processes Z(t) and Zjy(t) will always depend on q and on the 
{gk}i although we omit these dependencies in the notations. 

We want to prove that finite-dimensional distributions of processes Z(t) are 
limits of the distributions of processes Z^(t). 

More formally, introduce embeddings: 

l n : GT N <^Af, (Ai > • • • > Aw) i-^- (Ajv < • • • < Ai < Ai < Ai . . .). 

We also use the same notations for the induced maps A4(GTn) — > Ai(Af). Note 
that these maps are isometric in total variation norm. 
Cylindrical subsets of AT have the form 

U = {(>! < v 2 < . . .) €Af | v x G H u . . .,i/ fc G Hk) 

for aribitrary subsets Hi, ... , Hk of Z. 

The following statement is the main result of this section. 

Theorem 5.1. For every k > 1, the joint distribution of random vari- 
ables (in(Zn(0)), . . . , iAr(Zjv(A:))) weakly converges to the joint distribution of 
(Z(0),...,Z(fc)) as N ~> oo. 

Equivalently, ifAo,...,Af. are arbitrary cylindrical subsets of Af , then 

lim Prob{i A r(Z A r(0)) G A , ... , L N {Z N {k)) G A k } 

N— ¥00 

= Prob{Z(0)eA ,...,Z(k)eA k }. 

We start the proof with two lemmas. 

Lemma 5.2. Let /j be a finite measure on N ' , and let A be any cylindrical subset 
of N '. We have 

fi(A) = lim t N (E q {fi) N )(A). 

N— >-oo 

Proof. It suffices to prove the lemma for cylindrical sets of the form 

A = A b(1)! ... m ={v 1 <v 2 <---\v 1 = 6(1), ...,u t = b{£); b(j) G Z}. 

First, suppose that /i = S v for a certain v £ Af, then /i(A) = 1 if i/j = 
6(l),...,z^ = 6(£), and ^t(A) = otherwise. The statement of Lemma 5.2 
follows from Theorem 14.61 
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For an arbitrary measure we have 



fi(A) = [ 5"(A)»(dv) = f lim L N {E q {8 v ) N )(A)ii{dv) 

L N {E q (S v ) N )(A)fi(du) = lim L N (E q {fx) N )(A), 



(*) ,. 
= hm 

TV— >oo J \j * ^ ' ^ " y ' N — yoo 

where the equality (*) follows from the dominated convergence theorem 
Let us denote by I a the indicator function of set A: 



□ 



Ia(x) 



1, X £ A, 

0, otherwise. 



If /j is a measure, then l A [i stands for the measure given by 

(I A ti)(B)=ti(Ar\B). 

Lemma 5.3. Let \x be a probability measure on Af and let A be any cylindrical 
set. Then the total variation distance between measures L^(E q (filA)N) an d 
lAtN(E q ([j,)N) tends to zero as N — > oo. 

Proof. First, suppose that /x = 8 U for a certain sequence v £ AT. If v G A, then 
E q (lAH)N = E q (n) N , consequently, 



| l n (E q (I a^) n)-Ialn (E q {p) n) 



\{l-I A )L N (E q (ji) N )\\ = L N (E q (n) N )(A), 



where A — Af \ A. The right-hand side tends to zero by Lemma 5.2 
If v does not belong to A. then {I At 1 ) n 



\l N (E q (I A p,) N ) ~ I A L N (E q (fJ,) N )\ 



and 

I A L N (E q {n) N )\\ = i N (E q (p) N )(A) 
by Lemma [5. 2| 

To prove the claim for a general measure /x we observe the following property 
of the total variation norm. Suppose that W and V are measurable spaces, and 
/jv is a sequence of measurable maps from W to A4(V), such that ||/jv(iu)|| — 1 
and ||/jv(w)|| —> for any w € W. Then for any probability measure n on W, 
we have 



f N (w)ir(dw) 



w 



0. 



We obtain 

\i N (E q (I A pL) N ) - I A L N (E q (fj) N )\\ 

L N (E q (I A S") N ) - I A i N {E q {5 v ) N ) ),,(d,>) 



-> 0. 



□ 

Proof of Theorem \5.1\ Let us write P^ for the linear operator on AA(Af) corre- 
sponding to the kernel Poo(u, A] g(x)): 



(P^)(A) = / P x {u,A-,g{x))*{du). 
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Let P^r be the operator acting on measures of the form LN(A4(GTiy)) via 
the kernel P^(g(x)): 

AgGTjv 

Note that these operators are contractions in total variation norm. 
We have 

Prob{Z(0) g Ao, . , Z(k) G A k } = (I Afc P£-i • . • (/^P^CU^o)) • • 0(^0 
and 

Prob{L N (Z N {0)) G A , . . . , tiv(£iv(A0) G A fe } 

= (7 Afc P^- 1 ...(/ Al P^(/^ ))...)W- 
Applying Lemmas |5.2| |5.3| and definitions we obtain 



(I Ak P^...(I Al P^(I Ao Z ))...)(M) 

= lim tJV (P 9 (/ Afc P^- . . . (I^P^J^Zo)) . • . )n)(A0 

TV— »-oo 



= lim (7 Afc4jv (£;,(P»^ . . . (J^P* (I Ao Z )) . . . ) N ))(Af) 

iv— »oo 

= lim {I Ak P a N ^I Ah _,L N {E q {P%-*{. . . (I Al P^(I Ao Z Q )) . . . ) N ))(Af) 

N — i-rxn 



= • • • = lim (Ia^- 1 . . . {I M P% (I Ao Z Q )) . . .)(A0 

iv— >oo 

as desired. □ 
5.2 Correlation functions 

Let be a stochastic process taking values in subsets of Z (= point configu- 
rations in Z). For n > 1, the nth correlation function p n of is the following 
function of n distinct pairs (xi, 

/0 n (^i,*i;a;2,*2; ■ ■ • ;x n ,t n ) = Prob{a;i G Q(*i), . . . ,x n £ Q(t„)}- 

Recall that we identify GUV with iV-point configurations in Z. and A/" is 
identified with infinite subsets of Z that do not have — oo as their limit point. 



The following statement is a corollary of Theorem 5.1 



Theorem 5.4. The correlation functions of processes Zjsr(t) pointwise converge 
as N — > oo to the correlation functions of processes Z(t). 

Proof. Let us proof this statement for the first correlation function, for all other 
correlation functions the proof is analogous. Let pf (X, t) denote the first corre- 
lation function of the point configuration corresponding to Z^(t) and let p\(x, t) 
denote the first correlation function of the point configuration corresponding to 
Z(t). 

Choose e > 0. Let m be a number such that Prob{Z(£)i > m} > 1 — e. 
Since the distribution of i^{Z^{t))i converges to that of Z(t)\, for sufficiently 
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large N we have Prob{((,,/v.Zjv(i))i > m} > 1 — 2e. Now let k > x — m. Note 
that (ijv-2jv(i))i > m implies {bj^Zjs[{t))k+i + (k + 1) — 1 > x. Thus, 



Prob 



$y(i N Z N (t))i = x or ... or (i N Z N (t)) k + k - 1 = x\-px(x, t) 



< 2e 



and similarly for Z(t). 

It follows from Theorem |5.1| that for sufficiently large TV 



< 



< 



Pvob^(i N Z N (t))i — x or ... or (L N Z N (t)) k + k — 1 = xj 

— Prob|z(i)i =ior ... or Z(t)k + k — 1 = x\ 

Therefore, for sufficiently large N we have 
Pi(x,t) - pi(x,t)\ 

Pi (x,t) - Prob^(i N Z N (t))x = x or ... or (L N Z N (t)) k + k - lj 

Prob|(iArZAr(f))! = a; or ... or {l^ ZN(t))k + k — 1 = a; j 

— Prob|z(£)i = x or ... or Z{t) k + k — 1 = x j 

+ Prob|z(i)! = x or ... or Z{t) k + fc — 1 = xj — pi(x, t) 
< 2e + e + e = As. 
Since e is arbitrary, the proof is complete. 



□ 



Now we are in position to actually compute the correlation functions of 
Z(t). From now on we confine ourselves to the case Z = 6°. In other words, 
Z(0) = = (0 < < < . . .). As shown in [T71 Theorem 1.1], this implies that 
for every N > 1, Z N (0) = N = (0 < < • • • < 0). 

Proposition 5.5. For any n,N > 1, the nth correlation function p^[ of the 
process Z^{t) started from Zn(0) — Ojv admits the following determinantal 
formula: 

Pn(xi,h;x 2 ,t 2 ;...;x n ,t n ) = det [K N (x i ,t i ;x j ,tj)], 

i,j=l,...,n 



where 

K N (xi,t 1 ;x 2 ,t 2 ) 
1 



27Ti l^k+T fi/ iWltl 



>*2 



dw 



N-1 



-, (5.1) 



the contours C and C' N are closed and positively oriented; C includes only the 
pole and C' N includes only the poles q~ % , i = 0, . . . , N — 1 of the integrand. 
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Proof. See Theorem 2.25, Corollary 2.26, Remark 2.27 in [5], see also Proposi- 
tion 3.4 in To match the notations, one needs to set ag of [5] to be <? 1_£ , 
I > 1, set the symbol of the Toeplitz matrix F t (z) of [5] to gt(z _1 ), change the 
integration variables via £ i-> and shift the particles of |5] to the right by 
N. □ 

In what follows we use the standard notation 

oo 

O;g)co = ]]_{!- wq l ). 

Theorem 5.6. For any n > 1, £/ie ni/i correlation function p n of process Z(t) 
started from Z(0) = ftas i/ie /orm 

p n {xi,ti;x 2 ,t 2 ;...;x ni t n )= dot [AT(x i; ^; Xj, t,)], 

ij = l,...,n 

w/iere 

dz ffi! gt(w)(w;g) - 1 , 

w — z 



C is positively oriented and includes only the pole of the integrand; C goes 
from +ioo to —ico between C and point 1. 



Figure 1: Contours of integration for Theorem 5.6 C in blue and C in red 



Proof. It is sufficient to prove that for every (xi, ti; X2, £2), K {x\, t%; x%, £2) 
if (sci, ii; X2, £2) as AT — > 00. 



Denote by f^(z,w) the integrand in the second (double) integral in (5.1l. 
Note that if N is large enough, then for every w, \z 2 /n(z,w)\ goes to zero as 
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\z\ — > oo. Consequently, the replacement of the contour of integration C' N by C 
does not change the integral. 
Observe that the integral 



dzf N (z, w) 

a 

converges uniformly in w G C because of the rapid decay of f^{z,w) when 
z — > ±ioo. Moreover, the functions fpf(z,w) uniformly converge on C x C' as 
N — > oo. Therefore, 



Ilt'o 1 9t(w) ^(w) z** 1 



dw f dzf N (z,w)~^ <b dw r 
c Jc> Jc Je> Ut = o 9t{z) <t>( z ) w w - z 



□ 



5.3 First coordinate of the process 

In this section we give an independent interpretation for the evolution of the 
first coordinate of Z(t). Similar interpretations also exist for the evolutions of 
first k coordinates for every k. Theorems of this section are based on the results 
of [5, Sections 2.6-2.7] and we are not giving their proofs. 

Although all constructions make sense for general processes introduced in 
the previous sections, for simplicity of the exposition we restrict ourselves to 
homogeneous Markov processes started from the delta measure at 0. 

Denote by X g {t) a discrete time homogeneous Markov process on Af with 
transition probabilities P^iu, A; g(x)) started from the delta measure at 0. 

Also let 3^ 7+ . 7 _(t) be continuous time homogeneous Markov process on Af 
with transition probabilities (u, A; exp(£(7 + a; + 7_/x))) started from the 
delta measure at 0. 

Denote by X g {t)\ and y i+ , 1 _{t)\ projections of these processes to the first 
coordinate (i.e. we consider the position of the leftmost particle). Note that 
X g {t)i and 3^7 + ,7_ (t)i do not have to be Markov processes. 

Similarly, denote by X^ g {t)^ and 3^/v !7+j7 _ (t)jv the projections of processes 
■%N,g(t) and 3^Af.7 + ,7_ (t) to the iVth coordinate (which again corresponds to the 
position of the leftmost particle). 

As in [5], we introduce a state space S of interlacing variables 

s = (kt .. c | 4% < x™- 1 < *f } ■ 

^ m— 1,2,... J 

We interpret x" 1 < x™ < ■ ■ ■ < x™ as positions of m particles at horizontal line 
y = m. An example is shown in Figure [2] 
Introduce projection maps 

n N : S — > Mat = GT^r, 

{Xk)rn>l, l<k<m ^X? + N-1<X?+N-K'"<X% + N-1. 

Let % be a set of decreasing sequences of integers: 

H = {yi > 2/2 > 2/3 > ••• I y l G Z}. 
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Figure 2: Interlacing particles 



Denote by II the projection from S to H: 

n(K 1 }) = x\>x\>x\> .... 

Finally, let TV* denote the map from S to ZU {— 00} given by 
n°° (W}) = Jim x N + N —1= lim ^({xj 1 })^], 

where by ( ■ )[1] we mean a coordinate of the leftmost particle (it corresponds 
to iVth coordinate in GT/v)- 

An algebraic formalism which leads to a family of Markov processes on S was 
introduced in [5J Chapter 2]. Among the processes in [5] there were processes 
X g (t) and D^ 7+)7 _(t) such that the projections irff(X g (t)) and 7rjv(3^ + i7 _ (i)) 
coincide with Ajv, s (i) and 3^v,7 + , 7 _ respectively. Moreover, the projections 
II(A^(t)) and n(3^ 7 7 _(i)) are also Markov chains that we explicitly describe 
below. 

The following theorem explains the connection to our processes. 

Theorem 5.7. Finite dimensional distributions of the stochastic processes 
U°°(X^(t)) andU°°(y^ )7 _(*)) coincide with those of X g {t)x and y i+n _{t) 1 . 

In other words, Theorem |5 . 7| states that the first coordinate in the stochastic 
process X g (t) (or y i+tl _{t)) evolves as the limiting value of coordinates of the 
Markov process Il(xf{t)) (or II(y 7+l7 _ (t))) . 



Note that, in particular, Theorem 5.7 guarantees that if {x^} is distributed 
as Afjy,g(t)jv or 3^Af,7 + , 7 _ (4)jV) then the limit limjv-yooOsi — N + 1) is almost 
surely finite. 

Proof of Theorem \5. 7| The proofs for and 3^ 7+ , 7 _ (t)i are the same, and 

we will work with X g (t)i. 

Theorem |5.1| implies that finite dimensional distributions of the processes 
XN,g{t)N converge as N — > 00 to finite dimensional distributions of X g (t)\. 

Therefore, 

Probjn 00 ^/^)) g4- ■ ,n°°(A' g 5 (t fc )) e A fc } 

= Prob{ lim n N ({xT(ti)})[l] €A U .. . ,71* ({*?(**)})[!] € A k } 
= lim Piob{X Nj g(ti) N e Ai 9 . . . , X Ng (tk)N € A k } 

N—t-oo 

= Prob{^ g (ti)i 64..., # fl (f*)i G 

□ 
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It remains to describe the processes U(Xf{t)) and H(3^ +i7 _(t)). We have 
X£(0) = y!j +n _ (0) = (0 > -1 > -2 > -3 > . . . )• We view coordinates of the 
process as positions of particles in Z. Thus, Xg(0) = 7 _ (0) is the densely 
packed configuration of particles in Z<o- 

Let us start from H(Xg(t)). The description depends on g(x), and we discuss 
only the cases g(x) = 1 + (3x and g(x) = 1 + j3x~ ; for other possibilities we 
refer to [FJ Section 2.6]. 

Given a configuration {yi} at time moment t, to construct a configuration 
at moment t + 1 that we denote by {z{\, we perform a sequential update from 
right to left. 

For g(x) — l+f3x. the first particle jumps to the right by one with probability 
Pi := [3/(1 + f3) and stays put with probability 1 — p\. In other words, z\ = y\ 
with probability 1 — p\ and z\ = y\ + 1 with probability p\. For particle 
number k we do the following: If yk = z^-\ — 1, then this particle is blocked 
and we set Zk = yk- Otherwise, the particle jumps to the right with probability 
Pk := q 1 ~ k (3/ (1 + q 1 ^ k (3) and stays put with probability 1 — p^. 

For g(x) = l+/3/x, the first particle jumps to the left by one with probability 
pi := [3/(1 + f3) and stays put with probability 1 — p\. For particle number k 
we do the following: If y^ = Zk-i, then this particle is forced to jump to the left 
and we set = yk — 1 (one might say that particle number k was pushed by 
particle number k — 1). Otherwise, the particle jumps to the left by one with 
probability pk — q k ~ 1 j3/(\ + q k ~ 1 (3) and stays put with probability 1 — p^. 

Observe that the above two update rules are closely related. Indeed, l+/3x = 
(f3x)(l + From the probabilistic viewpoint, this equality means that 

one stochastic step for g(x) = 1 + f3x is equivalent to the composition of the 
stochastic step with g(x) = I + f3~ x /x and deterministic shift yi — > yi + 1 for 
all i > 1. 

For continuous time Markov processes 11(3^ 7 _ (t)) the interpretation is 
quite similar. Particle number k has two exponential clocks, "right clock" and 
"left clock" with parameters q 1 ~ k j + and q k ~ 1 j^, respectively. These two num- 
bers are intensities of the right and left jumps. When the right clock of a particle 
number k rings, it checks whether the position to its right is occupied (i.e. if 
yk = yk-i — If the answer is "Yes" then nothing happens, otherwise the 
particle jumps to the right by one. When the left clock of the particle number k 
rings, then this particle jumps to the left and pushes the (maybe empty) block 
of particles sitting next to it. 

The processes above with one-sided jumps are versions of the totally asym- 
metric simple exclusion process (known as TASEP) and long range TASEP, cf. 
|42j . [50] |31j . The processes with two-sided jumps was defined and studied in 
g] under the name of PushASEP. 

6 The Feller property 

In this section we show that transition probabilities 

Poo(m; A; exp(t(~f + x + J-/x))) 

of the Markov process y i+ n _ (t) satisfy the Feller property. The notion of a 
Feller process makes sense only for the processes in a locally compact space, 
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and this is the reason why we embed AT into the bigger locally compact space 
Af as in Section above. 



6.1 Extended boundary 

The aim of this section is to identify the local compactification Af of TV with a 

boundary of the sequence of sets GTat and links P?f(X — > /x) with the sequence 
of parameters ^ = q 1 ^ 1 , i > 1. These sets and links were introduced in the 
second part of Section [3] 

Let M = (M , Mi, . . . ) be an element of hm M{GT N ). We say that M is of 

(k) 

class fc if for N > k the support of Mm is a subset of GTjy'. We say that M is 
of class oo if for any N the support of Mjv is a subset of GT^' . 

Proposition 6.1. For any M E ]^mA4(&T]y) there exist unique 
M^M^M 1 ^/ 2 , • • • e ]pnM(GT N ) such that AP is of class i and M = 
M°° + M° + M 1 +M 2 + .... 

If M is a nonnegative then M°°, M°, M 1 , M 2 , . . . are also nonnegative. 

Proof. Decompose Mjv into the sum 

M N =M^+M^ + ... + M^ 

with 

supp(M<*>) c GT$. 
Clearly, such decomposition is unique and 

||Mv|H]T||mW||. (6.1) 

k 

Note that 

N N — M N-1 

for k < N - 2. Set M% = M$ for N > k, M% = p£ +1 M£ +1 , M^ 1 = T{m^ 
and so on. One proves that for every k, M k = (M^)„>i is of class k. 

Furthermore, ( |0 yields that J2l=o \\ m n\\ ^ ll M r+ill < ||M||. Conse- 
quently, the sum XX M fc is well defined. Set M°° = M - J2T=o Mk - Tt 
follows that M°° is of class oo and. thus, 



M = M°° +^M fe . 

fc=0 

The uniqueness of the decomposition and nonnegativity are immediate. □ 

Now we fix k and aim to describe the set of all class k elements of 
hjn A4(GTn). Observe that these are just elements of ]^mA4(TN) with rV — 
( grjp(mm(fc,Af)) w ^ links given by the restrictions of matrices P^i^ ~* /■<)■ Thus, 
we may use Theorem |4.2| Therefore, it remains to identify the set of all extreme 
coherent systems, i.e. the extreme points of lim jM„(GT^" n ^ ,Ar ^). 
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Theorem 6.2. The extreme points of M p (&T^ ln ^ k,N ^ ) are enumerated by 
signatures A € GT^. Let £ x be an element o/h^n A / f p (GT^ m "^ fc ' jV ' > ' 1 ) correspond- 
ing to A. ITien /or N > k, the Schur generating function of the measure £ N 
is 

^ ^'' vttrC ^ n G * 



n • 

71=^ + 1 



and /or Af G {1, . . . , k — 1} we have 
S(xi,...,x N ;£ N ) 



where 



n g^.^w-v..,^) 



For k = 0, £/ie unique element of ^im A / f p (GT^ rLm( ' fc ' jV ^) is £ , smc/i that for 
every N , Efj is the unique probability measure on the singleton GT^ . 

Proof. The case k = is obvious. We use the general ergodic method to prove 
this theorem for k > 1 (see [HI], |36| Section 6] and also |12[ Theorem 1.1]). 
Choose A € GT^ . There exists a unique system of measures M-j\ . . . , M N 

such that for every i < N, Af. ( A is a measure on QT<- min(k ' N)) ] = ~p\ +1 M^ +1 , 
and M N is the delta-measure at A. We call such a system of measures a primitive 
system of A. The ergodic method states that every extreme coherent system 
M G Urn J Mp(GT^" m( ' fc ' jV ^) is a weak limit of primitive systems. In other words, 
there exist sequences Ni — > 00 and A' G (j r [p^ r nm ( A '>'' c )) sucn that for every N and 
every /1 G GT^ 1111 ^''^, we have 

Mjv(/i) = lim M N \fj,). 

i~ > oo 

According to the definition of links p\ + i-, the Schur generating function of 
the measure M N has the following form for Ni > N > k: 

saH^S '■•-9 ) n=N+1 

Proposition 4.1 of [17J implies that the weak convergence of M N (/i) as i — > 00 
is equivalent to the uniform convergence of the Schur generating functions of 
N to those of Mm ■ 



measures M M to those of Mm . Observe that 



JVi 00 

n n °" 

n=JV+l n =N+l 

uniformly on as i — ¥ 00. Therefore, functions 

S X i(x!, ...,X k ) 
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should also uniformly converge. But this happens if and only if the sequence 
of signatures A 1 stabilize to a certain A <E GT^. Then measures Mjy converge 
precisely to £jj. 

Thus, coherent systems £ x contain all extreme points of 
lim A / f p (GT^ nin ^ fc ' Ar ' ) ' > ). It remains to prove that all of them are indeed 
extreme. But this follows from linear independence of the measures £jy 
which, in turn, follows from linear independence of Schur polynomials 
s\(xi,...,x k ). □ 



Now we are ready to describe the map E from VW(TV) to limM(GTjv), cf. 



Definition 4.1 Let 7r be a finite signed measure on TV. There is a unique 



decomposition 



7T° + 7T 1 + 7T 2 



such that supp(7r fc ) C GUV C TV and supp(7r°°) C TV C TV. 

By Theorems |4.2| and |4.6[ ir°° corresponds to a unique element of 



hmTW(GTTv) which can be viewed as an element M°° of limA4(GTj^). Sim- 
ilarly, by Theorems 4.2 and 6.2 7r fc corresponds to a unique element of 
^M(GT^ in(k ' N) ) which can be viewed as an element M k of hm .A/t(GTjv). 
Note that 

||M°°|| + ||M°|| + \\M l \\ + ••• = \y°\\ + || || + Ik 1 !! + ■ • • < oo. 
Therefore, we may define 

M := M°° + M° + M 1 + . . . . 

Set E(w) = AI. 

Theorem 6.3. The set TV and the map E satisfy the first 3 conditions of 



Definition J^.l 



Proof. 1. Proposition 6.1 implies that E is surjective. E is injective by the 
construction. 

2. E is a direct sum of norm 1 operators, thus it is a norm 1 operator. 
Furthermore, E is a bijection, thus E~ x is also bounded. 



3. Second part of Proposition 6.1 guarantees that E maps A4 p (Af) bijectively 
onto ^im M v (GTjy), 

□ 



As for the fourth condition of Definition |4.1| its proof is nontrivial and we 
present it as a separate theorem in the next section. 

6.2 The topology of the extended boundary 

In this section we prove the following theorem. 
Theorem 6.4. The map 

v -> £ v = E(S V ) 

from TV to hmTW(GTjy) (equipped with weak topology) is a homeomorphism on 
its image^ 



Note that the image of this map consists of the extreme points of lim Alp(GTjv)- 
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We present a proof in a series of lemmas. 

Let a[1..7V] be an element of Mat, i.e. a[1..7V] is the TV-point subset of Z 
consisting of points a[l] < a[2] < ■ ■ ■ < a[N]. Set 

A%( Xl ,...,x N )=Alt(x^...x^)= ]T (-irx?°W...xf^. 

aeS(N) 

Lemma 6.5. The following factorization property holds: 

s x ( Xl , . . .,x N ) = sx N _ k+1 ,...,x N (xi,..., Xfc)s( Al+fc! ,.. ;Ajv _ fc+fc )(x fc+ i, . . . , Xn) 
s x (l,...,q 1 - N ) " Ui=i,...,kUj=k+i,...M x i ~ x i) 
x 11, : ,11, v'V ' T ' = 

s\ N _ k+1 ,...,\ N (l, • • • , q 1 - k )s(\ 1+ k,...,\ N - k +k)(<i~ k , ■■, q 1 ~ N ) 

where if we expand Q as a sum of monomials 

Q= E 



+ Q, (6.2) 



c 



mi,...,m N mi 



J N ! 



mi r ..,mjv 



^ | croi ,...,m W | 1 r ni ... g (l-JVW <R ( k ,\ N _ k ,\ N _ k+1 ), 
mi,...,mjv 

and /or any fixed k and bounded Ajv-fe+i, -R — > as Ajy-fe — > oo. 

Remark 1. The statement should hold for more general sequences of £'s, 
but our proof works only for geometric progressions. 

Proof. We have (the first equality is the definition of the Schur polynomial) 

AH(t Xn t Xn - 1+1 ■ ■ ■ r Al+Ar ~ ls | 

(-l) N ( N -V/ 2 Sx ( Xl ,...,x N ) = — [Xl * 2 



Ui<j(Xi-Xj) 

Alt{x\ N ■ ■ ■ xl N - h+1+k - 1 )Mt{xl'll h+k ■ ■ ■ x x N 1+N - 1 ) 

Yii<j( X i ~ X j) 

Let us estimate the remainder B. By the definitions 

\ H1) +N-j(l) X Hk) +N-j(k) 



+ B. 



E 



{j(l),...,m}^{N-k+l,...,N} Ui<j<k( x i x o) 



X 



(X+5)\{X j(1) +N-j(l),...,X j(k) +N-j(k)}, . 
A N-k {Xk+1, ■ ■ ■ ,X N ) 

(6.3) 



rifc<i<j( a; i x j) 

i 



rii=i,...,fe ri i j=fe+i,...,jv( a '« x j) 

where (X + S)\ {A j(1) + iV - j(l), . . . , A j(fe) + TV - j(fc)} stands for the (N - k)- 
element subset of Z which is the set-theoretical difference of A + 5 = {Ai + TV — 
1, A 2 + N - 2 + . . . , A w _! + 1, A w } and {A j(1) +N — \ j{k) +N- j(k)}. 
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Denote the three factors in (6.3 1 by B\, B 2 and -B3, respectively. For any 
function L on T/v that has a Laurent expansion, let Est(L) be the following 
number: Decompose L into the sum of monomials 



m 

and set 



jrai,...,m N mi ~™»JV 



• • • X 



N 



mi,...,m;v 



Est(L) = |£ mi '-' mjv |l mi • ■ ■ 



mi ,....rn,N 

Clearly Est{LiL 2 ) < Est^Est^). Therefore, 



E'st(B) < Est(B 1 )Est(B 2 )Est(B 3 ). 

In what follows constk denotes various constants depending solely on k. We 
have 

k 

Est(B x ) < const k Y[ g( 1 -™)(^(m)+^-^(™)). 

m— 1 

Observe that i?2 is a Schur polynomial. As follows from the combinatorial 
formula(see e.g. Section 1.5 of (32] ) , these functions are sums of monomials with 
non-negative coefficients, thus Est(B 2 ) = B 2 (l,q~ 1 , ■ ■ ■ ,q 1 ~ JV ). 

Finally, decomposing the denominators in B 3 into geometric series and then 
converting them back, we get 

k \ N 

Est(B 3 ) < ( J] n (q 1 -! - q 1 - 1 )- 1 < const k J[ q^ k . 

,i=lj=k+l,...,N J j = k+l 

Now set a[N -j + l] = \ - + N- j. We have 



N k 



Est(B)< const k \ Yl q^-^Ylq^- 1 ^ 

{j(l),...i(fc)}/{l,...,fc} \j=k+l i=i 



.a[l..N]\{a\j(l)],...,a\j(k)]} , _fc 1 _JVx 



1 



.U[±..JVj u u L j WJ ,...,u[j WJJ - / 

(9 ""^ 'u l<3 < N -k(i- k - i+1 -i- k - ]+1 ) 

Let C(j(l), . . . , j(k)) denote the right-hand side of the above inequality. 
Then 

cm),...,m) A g(1 -o W (o]- W ) ^-f ] ^ b(1)1 -- ab - (fc)1 >( g -*, , g 1 -") 
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For any increasing sequence b[l..N — k] we have 

Af-?~ k] (<r fe ' • ■ • . I 1 '") = q (i-")m+---+b[N- k] ) A b[i : .N-k] (1) _ _ ) qN - k - 1} 

fe) q (l-N)(b[l]+-+b[N-k]) TT ( g &[»] _ g 6L/]) 
i<j<(N-k) 

= „( 1 - Ar )( h [ 1 H h6[A r -fe]) + (W-fe-l)b[l] + (Ar-fe-2)6[2]H |-06[jV-/s] 



x J] (i-gWMM) 



i<j<(N-k) 

= q (-k)b[l] + (-k-l)b[2] + - + (l-N)b[N-k] TT (1 _ g^bl-fcW)^ 

i<j<(N-k) 

where the equality (*) is the Vandermonde determinant evaluation. 
To analyze A f^]\{aU(i)],---A3(k)]} we think of the get 

b[l..N-k] :=a[l..N]\{a\j(l)},...,a[j(k)}} 

as of a small modification of the set a[k + 1..N]. Note that under this 
modification only finite number of members (up to k) of the sequence b[i] 
change. Using the finiteness of rin>i(l — 9™)' one easily sees that the 'mod- 
ified' product rii<i<j<(Ar-fe) (1 — q b bl~ fc M) differs from the 'unmodified' one 
Uk+i<i<j<N(^ ~ g°b']- a [*]) by a constant that is bounded away from and 
oo (note that k is fixed, while N can be arbitrarily large). 
Hence, 

A a[l..N]\{a[j(l)],...,a[j(k)]}( n -k „1-N) 
A N-k \1 ■•■•^ I 

^(r" 

< COnst k ■ q (-mi] + (-k-l)b{2]+- + a-N)b[N-k] q ka[k+l]+- + (N-l)a{N] ^ ^ 

The next step is to estimate the exponent of q above by 
k(a[k + 1] - 6[1]) H + (N — l)(a[N] - b[N - k}) 

N-k 

>k^2(a[m + k}-b[m})>k ^ (a[j (m)} - a[k}) . 

m=l m:j(m)>k 

Hence, 

C ^C{l "^k)^ ~ const k f[q {1 - l){a[j{i)] - a ^\ k ^^(^> ki ^ 

< const(k, a[k})- ^*«»)>* o[j(m)1 . 

Here and below we use const(k, a[k]) to denote any constant depending only on 
k and a[k}. 

Thus, the sum over all . . . ,j(k)} ^ {l,...,fc} can be bounded by 

geometric series and 

Est(B) < C(l,...,k)-const(k,a[k})q a ^ k+1 l 



Hence, substituting the definition of C(l, . . . , k) and using (2.5), 
B 



Est 



N k 

(l-i)(A JV+ i_ < +»-l) 



<const{k,a[k])q a ^ JJ g (i " 1)fc JJ? 

j=jfe+l i=l 

y „(-fe)(Ajv-_ fc +fc)+-+(l-JV)(Ai+iV-l) TT 1 Z ff A '~'~ Aj+J 

y 11 „-fe-*+l _ a -k-j+l 



x g 



l<i<j<N-k y y 
(JV-l)(Ai+-+A JV )„-OA 1 (JV-1)A W TT I - q° 1 



n r 



_ — /rAj— i-Aj+j 
l<i<j<Ar ^ 



= const(fc,a[fc])q a[fc+11 JJ gW" 1 )* 
j=fc+l 

AT ■_ 

1 — q J 



xn«- (M) n 

»=1 l<i<j<(Af-fc) 



„-fc-i+l _ q-fe-j'+l 



X 



1 - q3~ l j-r 1 - ^ 



TT TT tZH TT 

11 11 1 _ oAj-i-Aj+j 11 ]_ _ oAj-i-Aj+j 

i=l j=jV-k+l y (AT-fc+l)<i<j<AT y 



iV W 

< CO r^(fc, a [fc])<f[ fc+i i n 9 (i " i)fe n«" ( ' _i)2 n « fe+ ^ 1 

j=fc+l i=l l<i<j<(Ar-fc) 

N N N 

= const(k,a[k])q a ^ J[ q^ k J[ q~^ 2 J[ gO" 1 )^-*) 
j=k+l i=l j=k+l 

= const(k,a[k})q a ^ k+1 } 
The statement of the lemma immediately follows. □ 

Lemma 6.6. Let P be a signed finite measure on GTjy with Schur generating 
function S: 

* £™$b^r 

The total variation norm of P can be estimated as 

||P|| < const N Est{S). 

Proof. As before, we use the notation const^ below for various positive con- 
stants depending on N. We have 

Est(sl[( Xi - Xj ))= J2 |P(A)I Est(Alt(x^^... X ^)). 

\ i<j J AeGTjy ^ ' ' ' ' ' ^ > 

Note that 

Est Ultix^- 1 • • • 4")) > q^ql-mN-i+i) . . . q (i-N)(X 1+ N-i) _ 
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Also, using (2.5 1 



S A (1, ■ • • , q~ N ) < 1 q 0*N q (-l)(^-i+l) . . . q (l-N)(X 1+ N-l) 
' ' ~~ const n 

Therefore, 

Est I S Y\i x i ~ x j) ^ constpf | J-*(A) j = const jv||P||. 

\ i<j J AGGTjv 

It remains to observe that 

Est Sj\[xi - Xj) < Est{S)Est Y[(x t - Xj ) < const N E st(S) . 
V i<3 I \*<J / 



□ 



Lemma 6.7. Le£ z/™ 6e a sequence of points ofJ\f C M converging to v € M\M . 
Then £ v weakly converges to £ v . 

Proof. Since v g Af\Af, there exists fc > 1 such that ^ € GT^ C A/". First, take 
m > k. Recall that the measure £" n is supported by GT m and S(x\, . . . , x m ; 
is its Schur generating function. 
Theorem 1.3 of [T7] yields 

S(xi, . . . , x m ; £^ ) = Jim S 1 ^ (xi, . . . , x m ), 

N-too 

where 

„„™ , v _ S^,...,^)^!) ■ ■ ■ ; x m, g m ,---,q ) 

^ 1*1,..., *m; - , -jz^ 

Since v n — > we have v^+i ~~ ^ °°- wm l e ff , . . . stabilize. Thus, we may 
use Lemma l6~5l 

so», ...,!/?) fai, • • ■ , im, q~ m , • ■ • , q 1_N ) 



•S 



S^n^n^Xx^.^Xk) Ui=l,...,kUj= m +l,..., N (^ ^ 

s 



( ^,..., y n)(i,...,g i fe ) n 4= i,..., fc n J=m +i,...,Ar(^-9 1 3 ') 

s(^+fe,...,i^ +1 +/£)( :z: fe+i I ■ ■ ■ , x m,q~~ m , ■ ■ ■ ,q 1 ~ N ) 

S(^+k,...,^ +1 +k)(q~ k , ■ ■ • , q^ N ) 



11, i „.V ' '/' * 

rii=i,...,fc rij=fc+i,...,m( :i: j _ x j) 



Q. (6.5) 



Note that the measure with Schur generating function S 1 ^ (x\, . . . ,x m ) is 
supported on signatures A such that \ m -k > v k+i- Indeed, this readily follows 



from formulas ( |3.1| and (3.2 1, see also Proposition 5.5 in |17j . 

Let us find the projection of the measure with Schur generating function 
(xi, . . . , x m ) to the last k coordinates of the signature (all other coordinates 



3G 



tend to infinity as n — > oo). We claim that the Schur generating function of this 
projection uniformly (in x's and in N) tends (as n — > oo) to 



SW,...,v»)(Xu...,Xk) II, i /.II 



j=m+l,...,N 



(<? 



l-i _ „l-j 



( v k< 



Iii=l,...,klij=m+1,...,N\ X * I 1 °J 



.6) 



To prove this claim, we first expand (6.6 1 into a finite sum of normalized 
Schur polynomials 

EN s x{ x l> ■ ■ ■ ! ^fe) , ^ 



SA(l,-.-,5 1 - fe ) 



It is clear, that for any e we may choose such a finite set A not depending on 
N, that Est(Qi) < e. 
We also expand 

S( w »+k,..., w jf +fc)(a*+i, ■ ■ • ,x m ,<T fc , . . .,g 1_JV ) 



3 (i^+fc,...,i^ +1 +/c)(<7 fe , • 



■ ,q 



l-N 



into the /wZZ sum of normalized Schur polynomials in variables x^+i, • ■ ■ , x n 
(with normalization in the point (q~ k , . . . , q 1 ^" 1 ) as 



E 



i s ll (q- k ,...,q 1 - m ) 



Observe that ^^u^ — 1. Note that, if n is large enough, then coordinates 
of all signatures in the support of u M has larger coordinates than those of A. 
Moreover, as n — > oo these coordinates tend to infinity. 



We substitute these two expansions into (|6.5| and use Lemma|6.5|yet again 
ion, for po 

E^E 



(in the reverse direction, for polynomials in m variables). This gives 

S\Ufj,{ x l, • ■ ■ i x m) 



l — m\ 



7) 



If we forget about Q2, then we arrive at a measure on GT TO which assigns 
to signature XL) fi the weight c\u^. Clearly, its projection to the lowest k 
coordinates assigns to signature A the weight c\, as needed. 

It remains to work out the impact of Qi. Let us list all the terms contributing 
to Q 2 . 



1. Q from (6.5 1. By Lemma 6.5 Est(Q) — > as n — > 00. 

2. Qi gives the term 



Qi = Qi 



i (v™+k,...,v% +1 +k)( Xk + 1 > ■ ■ ■ i X m,q 

s r ..n , u ..»+! , ,„)(q~ k , ■ ■ ■,q 1 ~ N ) 



■ q 



1-N\ 



>™+/c,..., t ^+ 1 + fc)^ 



11, : a 11, „.!'/' ' '/' * 

rii=l,...,fc Ylj=k+l,...,m( X i ~ X j) 



We have EstiQ'J < Est(Qi) < e. 
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3. When we use Lemma 6.5 the second time, we get a term Q A ' M for each 



pair (X,fj,). Thus, the total impact is 



A.// 



Since by Lemma 6.5 we have uniform bounds for Est(Q x,fl ) (here we use 
the fact that A is finite) and J2\ a l ? V c i V | < 1; thus, the contribution of 
these terms also tends to zero as n — > oo. 

Summing up, we have 

Est{Q 2 ) <s + R, 
where R — > as v n — > v (uniformly in N) . Consequently, 

Est(^Q 2 Yl ( Xl ~ X i)) - const rn(£ + R)- 
l<i<j<_m 

Now let H be a (signed) measure on GT m with Schur-generating function 
Q 2 . Using Lemma |6.6| we conclude that total variation norm of H can be 
bounded by 

< const m (e + R). 

Since e is arbitrary and R — > as n — > oo, the influence of H and, thus, the 



influence of Q 2 in (6.7) is negligible. 

We have proved that the Schur generating function of the projection to 
the last k coordinates of measure on GT m with Schur generating function 



(xi, . . . ,x m ), tends as n — ¥ oo to the function given by (6.6 1. Now sending 
N to infinity we see that the Schur generating function of the projection to 
the last k coordinates of measure £^ uniformly tends to the Schur generating 
function of £^. Next, we use Proposition 4.1 of [17: which yields that weak 
convergence of measures is equivalent to uniform convergence of their Schur 
generating functions. We conclude that the projection to the last k coordinates 
of measure £^ weakly tends to £^. On the other hand, all other coordinates of 
the signature distributed according to £^ tend to infinity as n — > oo. It follows 
that £^ — > £^. 

For m < fc, Lemma |6 . 5| yields 



„», S (l/ n_ v n ) ( Xl ,...,X m ,q m ,...,q 1 N ) 

"("wi-^Dv 1 ' ' ' ' q > 
_ Sjvn,...,^)^!, ...,x m , q~ m , . . . , g 1 "' ) Ui=i,..., m rij=fe+i,...,jv-(9 1_< - 



and 

lim S(x\, . . . ,x m ;£^ ) = lim lim S V N {xi, . . . ,x m ) = 

n— too n— too N^roo 

S( ^...., v?) (x u Xm , g-™, . . . , g^) uti n~ fc+ i(g 1 ~ i - t^l 

= S(x\, . . . , x m ', £ m )- 
By Proposition 4.1 of [T7] we conclude that £^ — > £^. □ 



-Q 
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Lemma 6.8. Let v n be a sequence of points of J\f\J\f converging to v £ Af\Af. 
Then £ v weakly converges to £ " . 

The proof is similar to that of Lemma |6.7| and we omit it. 



Lemma 6.9. Let v n be a sequence of points ofM\Af converging to v £ M. 
Then £ v weakly converges to £ v . 

Proof. Fix m > 1. Since v n — > v £ A/", for large enough n we have v n £ GTfc re 
with k n > m. Moreover, fc„ — > oo as n — > oo. Therefore, £^ is supported on 
GT m , and its Schur generating function is 

<j/ . c-w™N _ V^lr-'i^mig "V--;? 1 fc ") TT 

0^1,...,J, m ,C m ; — /j l-fc„) 11 9 ' 



Therefore, 



Jipi<S(xi, . . . ,x m ;££) = Jim ( ; 



S v n (X\ , . . . , X m) Q 71 ) 

= tS(xi, . . . , x m , £ rn ) 1 



where the last equality follows from Theorem 1.3 of [T7]. Using Proposition 4.1 
of [17] we conclude that £^ — > £^ n . □ 

Lemma 6.10. Let v n be a sequence of points of M C 77 converging to v £ J\f. 
Then £ v weakly converges to £" . 

Proof. See Proposition 5.16 of [17]. □ 

Lemma 6.11. Let v n be a sequence of points of 77 and v £ 77 '. weakly 
converges to £ u , then there exists to £ Z sitcft £/ia£ z/™ > to /or every n > 1. 

Proof. We start by proving that that there exists a constant c > such that 
£i ({^™}) > c f° r every n. > 1. More exactly, one can take c = (q; q)^. 

For v £ M this was proved in Lemma 5.15 of [17 . If v n £ GTfc, then 
the support of £\ n consists of numbers greater or equal then z/™. Therefore, 
x~ Vl S(x; £ i ) is a power series (without negative powers of a;) and £^ ({^™}) 
equal the value of this series at x = 0. 

We have 

/ —1 1 — k\ 00 



•=k+l 



Observe that s M+9 (xi, . . . , Xk) — (x±, . . . , Xk) q s fl (xi, . . . , x*;) and use this equal- 
ity for fi = is n , q = i/f. We obtain 

/ —1 1 — k\ 00 

^ SA (l,g-i,... )g i-fe) 11 ^ 
where \ — v n — is a signature with Afc = 0. Therefore 
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Using (2.51 we obtain 



sa(1,<Z -1 , 



,.,q 



n 



l _ g-Ai+Aj+i-j 



l-k\ 



i<i<j<k-i I — q* 

^ _ g-Ai+Aj+i-j 



,-|A| 



n 

l<i<j<fc 



1-9 



i=l 1 ~ 



-Ai+i-fe 



,i-fc 



fe-1 

Hi 



i-« 



fc-j 



fc-i 



i=l ~ q i=l 



>H0.-q h - i )>H(l-q i ). 



i=l 



Also 



Hence, 



fj GK0)-(l-g fc )(l-g fe+1 ) 



Now let C q {u) be a function on GTi = Z U {oo} that vanishes at all points 
> g and equals 1 at all points < q. Choose q so that 

C q {u)£»{u) < c/2. 

ueGTi 



Note that C q is a continuous function on GTi. Thus. 

Y, C q {u)£l{v)= lira £ C»£r> 

* ' 7?.— VOO * » 



Hence, for large enough n we have 



But if i/J 1 < g then 



^ C,(u)fj'"(u)>fi'"(K})>c. 



This contradiction proves that > q. 

Now we are ready to give a proof of Theorem |6.4| 



□ 



Proof of Theorem 6.4 Let {v n } n >i C A/" and {/ € A/". Our goal is to prove that 
lim, i _ i . 00 v n = v if and only if weakly converges to f asn-> oo. 

First, suppose that i/ n — > v. Without loss of generality we may assume that 
either for every n we have v n € Af, or for every n, v n € AT \ A/". Also either 
v £ Af or is e Af\Af. Thus, we have four cases and they are covered by Lemmas 

Now suppose that £ v weakly converges to £ v . By Lemma 6.11 there exists 
m € Z such that v n £ D m for every n > 1, where 



Dr. 



{i/ e A" : v x > m}. 
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Observe that the set D m is compact. Therefore, the sequence v n has a converg- 
ing subsequence v nh — » v' . But then £ v — > E v and, thus, v = v' . We see that 
{V"} is a sequence in a compact set such that all its converging subsequences 
converge to v. This implies v n — > v as n — > oo. □ 



6.3 Feller Markov processes on the boundary 

In this section we prove that for any admissible function g(x), P 00 (w,A;(;) is 
a Feller kernel on Af. Moreover, we show that a Markov process on M with 
semigroup of transition probabilities Poo(w, A; exp(t(7 + a;+7_/a;))) and arbitrary 
initial distribution is Feller. 

The section is organized as follows: First we prove in Proposition |6.12 that 
for any N > and any admissible g(x), P jv(A — > fi; g) is a Feller kernel on GTjy. 



As a corollary, we show in Theorem 6.15 that Poo(u, A; g) is a Feller kernel on J\f 



Finally, we prove (Proposition |6 . 16 and Theorem 6.171 that Markov processes 
with semigroups of transition probabilities Pn(h — > A; exp(i(7 + a; + 7_/x))) and 
Poo(m, A; exp(t('j + x + 7_/x))) are Feller. 

Proposition 6.12. For any admissible function g(x), i"V(/i ~^ A;g) is a Feller 
kernel on GT^r. 

First, we prove a technical lemma. 

Lemma 6.13. For any elementary admissible function g(x), the transition 
probabilities -P/v(/i — > X;g(x)) admit an exponential tail estimate: There exist 
positive constants a± and ai such that 



A; q(x)) < a\ exp(— a-i max |Aj 

" y Ki<N 



Mil)- 



Proof. Proposition 2.4 gives 
/ N 

n 



p N {ii 

where 



giq 1 - 1 ) 



dot 

i,j=l,...,JV 



sx(l, 



kez 



If g{x) = (1 + j3x ), then Pn((J- —> A;<7(x)) = as soon as |A$ — > 1 for 
any i, and we are done. 
For the remainii 



For the remaining two cases note that (2.5 1 implies 



It follows that 



< const N • gE^-^W-AO < const N ■ q~ N max * ^"^l 



,1-N\ 



det 

i,j=l,...,N 



C\,-; 



N 



<const N Q~ N maJE<|Ai_Wl n c ^-i-^ (i )+-(i)- 

creS(N) i=l 
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If g(x) = exp(7X ±:L ) then Ckr k — > for any r > as k — > oo. Note that 
maxi (A; — fj,i\ = m implies that for any permutation a there exists i such that 
|A* ~ Mo-(»)l > m - For this i the product 

— TV 2 maxi |A* — /i{L 

is exponentially small (in to). Therefore, each term 

N 

„-iV 2 roaxi |A<- W | TT 
i=l 

tends to zero exponentially fast as max^ |Aj — — > oo and we are done. 

Finally, if g(x) = (1 — ai -1 ) -1 then -P/v(/i — > ^i9( x )) — unless ^ > Xi for 
1 < i < N. In the latter case, 

— — < const. 

s Al (l,...,g 1 w ) 

Expanding the determinant 

JV 



det 

i,j=l,...,N 



by the same argument as above we see that each term in the sum tends to zero 
exponentially fast as max^ |Aj — Mil 00 ■ Therefore, 

' N -, \ r -i 1-JV 



\i=l 

tends to zero exponentially fast. □ 

Proof of Proposition^!^ Let h e C (GTV); we want to check that P^(g)(h) £ 
Co(GTjv)- By the definition of admissible functions, it suffices to check this 
property for elementary admissible functions g(x) — (1 + (3x ±1 ), g(x) — 
exp(7x ±1 ) and g(x) = (1 — ai -1 ) -1 . Moreover, since Co(GTjv) is closed and 
Pj!f(g) is a contraction, it is enough to check this property on a set of functions 
h whose linear span is dense. 

We choose the following system of functions (AG GUV, k>l): 



1, fi = X, 
0, otherwise; 



1) Hi > ■ ■ ■ , > Afe, /Ufe+i — Afc + i, . . . , f-tj\[ — Aat, 
0, otherwise. 



Let us start from h = a\(n). We have 

[~Pj(g)(h))(fi) = P^( P ^\;g) 



By definition, Pjv(jU —> X; g) = if fi € GTV \ GTV. Note that for a 
sequence /i™ of elements of GTV, A 4 " — > oo (in topology of GTV) means that 
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jiijy — > — oo, while /i" ->/i£ GIV \ GT^ implies /j,™ — > +00. Thus, to show 
that P^(g)(h) £ Co(GTjv) we should prove that if fi n is a sequence of elements 
of GTjv such that either //.jy — > —00 or /i™ — » +00 as n — » 00, then PN([J- n — > 
A; 5) = P/v(u™ — > A; g) — > as n — > 00. But if u'^ — > —00 then |/x^ — Aw| — > 00 
and we may use Lemma |6.13 If /1" — > +00, then — Ai| — > 00 and Lemma 



|6.13| also provides the required estimate. 

Now let h = b\ t k{p)- The fact that ^-F^rC?)^)) (") — > as /i — ^ 00 (in other 



words, as fijy —00) again follows from Lemma 6.13| However, we need an 
additional argument to prove that the function P^(g)(h) is continuous on GTjy. 

Let us prove that if [i n is a sequence of elements of GT^r converging to //, 
then (^P^(g)(h)j(n n ) — > (^P^(g){h)j{n)- One readily sees that there are two 

principal cases (the result for all other cases is a simple corollary of the results 
in one of these two cases) : 



• (S 1 £ GT N C GT N and fx £ GT N C GT N ; 

• /j, n £ GTV C GT^r and £ GTjy \ GT^r, i.e. //1 = \ii = ■ ■ ■ = « m = 00 
while [i m+ i < 00 for some 1 < to < N. 

In the former case the sequence u™ stabilizes, i.e. /1" = (i for large enough n. 
therefore, (^P^(g){h)j(fi n ) — ^P^(g)(h)j (") for large enough n. In the rest of 
the proof we concentrate on the latter case. From the definition of convergence 
in GTjv we conclude that — > +00 for i = 1, ... to, and for n > no we have 
/i™ = /i{ for i = m+ 1, . . . , N. Without loss of generality we may assume uq = 0. 
If to > k then by Lemma |6.13| 

(pj(g)(h))( l x n )^0=(p* r (g)(h))^). 
Now suppose that m = k. We have: 

(^G?)W)(M")= E PN{» n ^v,9), (6.8) 

where B\ ik = {116 GT N : b\. k (u) = 1}. 

Note that in /i™ the first k coordinates are large while the last N — k coor- 
dinates are fixed, and choose an integral sequence r" such that r" — » +00 and 
fi r k l — r n — > 00. Let 

A n = {v £ GT N I u x >r n ,...,v k > r n }. 
For every n we divide the set B\.k in the sum ( |6.8[ ) into two disjoint parts 

E ^(m"^^;p)= E E *W->";j7) 

ueB XM r/es^fc ^eBj_ fc 

where 

^A, fe = SA, fc n^l", BS ifc = B A)fc \^. 

Observe that, as follows from Lemma |6.13| 

E P N (^^^g)^0. 
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In the remaining sum 



P N ( t S l ^v;g), 



(6.9) 



we use Lemma |6.5| for every term, i.e. for 

1 



N 



\i=Ti g ( q > j ,J=1 '-' JV 



"ui—i—^+j 



As we have shown in the proof of Lemma |6.13| for all elementary admissible func- 
tions, the coefficients decay rapidly as k grows. Therefore, the determinant 



of the matrix 



'^T-i-Vj+j 



factorizes. We conclude that 



P N (fi n ^ v\g{x)) = (l + o(l)) 

/N-k j \ 

x TT / -s det 

V 0(0 ') / M=fe+l,...,JV 



'(Kfc+1,...,!/Jv) 



t N 

n 

Vi=iV-ft+l 



det 



C Vi -i-n™ +j 



S(i/i+JV-fe,...,i/ fc +JV-fe) w >•••># J 

s (j[1?+ iv-fc,... )/1 »+iv-fe) (<7 fc_Ar , ■ • • , <7 1_JV ) ' 



(6.10) 



where the term o(l) uniformly tends to zero as n — > oo. 

Note that the second line in (6.101 is the transition probability 



Pn[H^ (oo,...,oo,A fe+ i,...,Ajv);5(x)J = (P^(g)(h))(iJ,) 



As for the third line, observe that Lemma [2.2| implies 



A' 



E TT -tt^ det 

Oi,...,i/ fc )eGT fc \i=JV-fe+l MW V 



"ui-i-fi'j+j 



l„k-N ri l-N\ 
s (^i + A r -fe,...^ JV +A r -fc) W ,-■■,<? J 



s (/x"+jV-fc,...,/i£+Ar- 



-k)(« 



fc-JV nl-JV 



= 1. 



Arguing as in Lemma 6.13 we conclude that we may replace the summation 
set by B" fe , i.e. as n — > oo, 



E ( II g r q i 

veB" ,, \i=N-k+l 



) det 




/ i,j=l,...,k 





s (ui+N-k,...,u N +N- 



-k)(q' 



k-N 



^- N ) 



s(^+N-k,...^i+N-k)(q k Ar ,---,9 1 N ) 

Summing up, we proved that 

(Pj(g)(h))(fx n )= ^W l ^;s) + E *M** n ->";s) 

= (Pj(g)(h))(ri(l + o(l)) + o(l)^(Pj(g)(h))(ri. 
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It remains to consider the case m < k. The statement in this case essentially 
follows from the case m = k. Indeed, decompose i?A.fe into a disjoint union 

B\,k = [J Bg >m , 

where the union is taken over the set 9 C GTat consisting of all 9 such that 
9i = Xi for i > k, 9i > Xi for m < i < k and 9i — max(Ai, 9 m+ \) for i < m. Now 
choose large enough s and denote 

e s = en{ve gt^ | v x < s }. 

We write 

eee s eee\e s 

Note that the set s is finite. Therefore, using the already proven case m = k 
we conclude that 

E (^G?)(^,m))(M n ) "> E (^(fl)(^.m))(M). 

To finish the proof it remains to note that as follows from Lemma 6.13| 

E (Pj(9)(be,m))(n n ) 

uniformly (in n) tends to zero as s — > oo. □ 



Let -P^^jv be a contraction operator from S(GTjv) to B(N) given by: 

where £" is the extreme coherent system corresponding to the measure 6" on 
AT. 



Lemma 6.14. P^ N maps C {GT N ) to C (N). 

Proof. Again it suffices to check the lemma on a set of continuous functions 
whose linear span is dense. We choose the familiar system of functions 



1, fx = \, 
0, otherwise, 



&A,k(/f) 



1, [J-i > X\, . . . , [ik > Afe, jtXfc+i — Afe+i, . . . , fiN — Aat, 
0, otherwise, 



and their finite linear combinations. Theorem 6.4 implies that -P^o_>.jv ( a -0 an d 
■^tx)->iv(^A,)fc) are continuous functions on AT. It remains to check that they 
vanish at the infinity. If v n is a sequence of elements of A/" tending to infinity 
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(i.e. escaping from every compact set), then uf — > — oo. Thus, we should check 
that if Vi — > — oo then 

Assume the opposite. Then there exists a subsequence {ni} such that £ V N (A) > 
c > for n — rig. Let ifr = v 111 — v™ 1 . Then is a sequence of elements 

of a compact set. Hence, {tp £ } has a converging subsequence. Without loss 

of generality assume that already ip e is converging, ip^ — > ip. Since, £% is a 

probability measure on GTjv and — » Sfj, we must have £^ (A — z/™') — > 0. 

Now observe the following property of measures £ v which was proved in |17[ 
Proposition 5.12]. For e G Z and ^ g A/" let i/ — e be a sequence with coordinates 
(y — e)i — Vi — ej. In the same way for A € G1V set (A — e)i = Xi — e. Then we 
have£j^ e (A-e) =fJ^(A). 

We conclude that (A) = f^(A - v^ 1 ) 0. Contradiction. 

The argument for the functions b\^k is similar and we omit it. □ 

Theorem 6.15. For an admissible g(x), Poo(g) is a Feller kernel on Af. 



Proof. This is an application of Proposition 4.4 Indeed, GT^r and J\f are locally 



compact with countable bases by definition. The property 2. of Proposition 4.4 



is Lemma [6 . 14| and the property 3. is Proposition |6.12| □ 



Proposition 6.16. A Markov process on GTjy with semigroup of transition 
probabilities Pjv(A — > /i; exp(t(j + x + 7_/x))) and arbitrary initial distribution 
is Feller. 

Proof. The first property from the definition of a Feller process is contained in 
Proposition |6.12| As for the second property, it immediately follows from the 
fact that exp(£(7 + a; + j-/x))) —> 1 as t — > and definitions. □ 

Now Proposition |4.5| yields 

Theorem 6.17. A Markov process on 77 with semigroup of transition probabil- 
ities P oc (m, A; exp(t(7 + a; + 7_/ir))) and arbitrary initial distribution is Feller. 



7 PushASEP with particle- dependent jump rates 

Fix parameters £i,...,(jv > 0, a, b > 0, and assume that at least one the 
numbers a and b does not vanish. Consider N particles in 7L located at different 
sites and enumerated from left to right. The particle number n, 1 < n < N, has 
two exponential clocks - the "right clock" of rate a£ n and the "left clock" of rate 
b/Cn- When the right clock of particle number n rings, it checks whether the 
position to the right of it is empty. If yes, then the particle jumps to the right 
by 1, otherwise it stays put. When the left clock of particle number n rings, it 
jumps to the left by 1 and pushes the (maybe empty) block of particles sitting 
next to it. 

Note that if Ci = q^ N C2 = q 2 ~ N ', Ov-i = <r\ Cn = 1, and the 
process is started from the initial configuration 1 — N, 2 — N, 1,0 then 
this dynamics describes the evolution of N leftmost particles of the process 



n(3v, ~,_{t)) introduced in Section 5.3 
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Let Pt(x\, . . . , xn I yi . . . ,yj\r) denote the transition probabilities of the 
above process. The probabilities depend on fx, . . . , a, b, but we omit these 
dependencies from the notation. In this section we study the asymptotic be- 
havior of particles as time t goes to infinity. In other words, we are interested 
in the asymptotics of P t (x\, . . . , x^ \ yi . . . , yjsf) as t — > oo. 

One easily checks that if for some r £ M > o and 1 < k < N we have Cfc < r 
and £k+i > r, ... , Cn > T, then at large times the first k and the last N — k 
particles behave independently. Thus, it is enough to study the case £jv < 
miii(£i, . . . , (jv_i). Moreover, without loss of generality we may also assume 
that £jv = 1- Thus, it suffices to consider the situation when we have h < N 
indices ri\ < n 2 < ■ ■ ■ < — N such that 



and £k > 1 for k not belonging to {n\, . . . , n/j}. Set D = {ni\i=\ h' 

To state the result on asymptotic behavior we need to introduce a certain 
distribution from the random matrix theory first. Let M n be a random n x n 
Hermitian matrix from the Gaussian Unitary Ensemble, see e.g. |33| . |15| . [T] 
for the definition. We use the normalization for which diagonal matrix elements 
are real Gaussian random variables with variance 1. Denote the eigenvalues of 
M„ by A™ < A£ < • • • < A™. Let M k be the top left k x k submatrix of M n , 
denote the eigenvalues of by < • • • < A^; . Finally, denote by GUE™ the 
joint distribution of the smallest eigenvalues of matrices Mfc. In other words, it 
is the joint distribution of the vector (A™ < A™ -1 < • • • < A}). 
Denote by <3 n (z), n £ Z, the integral 



where the contour of integration C2 is shown in Figure p] Note that &o(z) is 






Figure 3: Contour of integration Ci. 



the density of the Gaussian distribution: 




exp(-z 2 /2). 
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More generally, for n > 1 we have: 



z 



where -ffn is the nth Hermite polynomial. 

Proposition 7.1. TTie probability distribution GUE n of the smallest eigenval- 
ues of the GU E principal submatrices has density 

p(Vl, ■■■,Vn) = det[<Sk-i(yi)]k,l=i,...,n, Vi < 2/2 < ••• < Vn- 

Proof. We start from the formulas for the correlation functions of all eigenvalues 
of GUE minors found in [2T], [37], [25], [23]. Given a n x n GUE matrix, one 
constructs a point process on N x K which has a point (to, y) if and only if y 
is an eigenvalue of to x m top left submatrix. Let p n be the nth correlation 
function of this point process. Observe that the interlacing of the eigenvalues 
of the nested submatrices guarantees that 

Pn(l,y n ;2,y„_i; . . ■;n,y 1 ) ) y n > y n -i > ■■ ■ > y lt 

is precisely the density of the distribution GUE™. Then Theorem 2 in [37] and 
Theorem 1.3 of [22 yield that the density of GUE™ is 

det [K(i,y n+1 _i;j,y n+1 _j)], 

% , j — 1 ... n 

where (we use a more convenient for us integral representation for the kernel, 
which can be found after the formula (6.17) in [22]; note that we use a different 
normalization for GUE) 



J2 s - r exp( ., 



l — n^Ci J C2 



X 

fc=0 



where C\ is anticlockwise oriented circle around the origin and C2 is the line 
$lv = 2 from — ioo to +ioo. Changing variables z — —\/2v and w — —\T2,u we 
arrive at 



(-1)' 



. ie xp(^ 



(2ttC 

v ' fe=0 



00 

V / dwexp(-(w - w 2 /2) w fe -''© s _ fc _ 1 (n), 



where C3 is clockwise oriented circle around origin. Note that when k > n > r 
the integral J c dw exp(—£w — w 2 /2)w k ~ r vanishes. Thus, we may assume that 
the sum is over k = 0, 1, . . . , n — 1. 

Hence, K(r, £; s, n) is the matrix clement of the product of two matrices and 

.det [K(i,y n+1 - i ;j,y n+1 - j )] = (-l) n det [A(r,k)} det [<8 s _ fe (y n+1 _ s )], 

i,j — l...n r,/c — l,...,n /c,s — l,...,n 

where 

A(r,k) = — / exp(y n+ i_ r w - w 2 / 2 
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Again note that A(r, k) = unless r > k. Thus, the matrix [A(r, k)] is triangular 
and 



det [A(r, k)] = TT A(k, k) 

r,/c— l,...,n 

fc=l 

with 

. . 1 f dw , 2 ;o\ 

A \k,k) = -— — exp(y n+ i_ r w - w /2) = -1. 
2m J c ^ w 



To finish the proof it remains to note that 



det [<S s -k(y n +i-s)} = det [<S k - s (y s )}. 

k.s—l,...,n fc,s=l,...,n 



□ 



Now we are ready to state the main theorem of this section. 

Theorem 7.2. Fix parameters ({Ci}" = i 7 a, b) as above and let y\ < ■ ■ ■ < y^ be 
arbitrary integers. Denote by (Xi(t), . . . ,X/v(i)) the random vector distributed 
as 

P t (xi, ...,xn\vi,-- -,Vn)- 
The joint distribution of 

f X ni (t) - (a - b)t X n2 (t) - (a-b)t X nh (t) - (a - b)t 
V yj{a + b)t ' y/(a + b)t y/(a + b)t ' 

X ni (t) - X ni _i(t), . . . ,X 2 (t) - Xx{t); 

X n2 (t) — X ri2 ^i(t), . . . ,X ni+2 {t) — X ni+ i(t); 

X n h {t) ~ X nh -i(t), . . . ,X nh l+2 (t) - X nh l + 1 (t) 

converges to 

GUE\ 

x Ge(C 1 1 _ 1 ) x ••• x Ge(Cr 1 ) 

xGe^x-xGefC^), 

where Ge(p) is a geometric distribution on {1, 2, . . . } with parameter p. 

Theorem |7.2| follows from Proposition |7.1| and Proposition |7.3[ 
Proposition 7.3. Set v = a — b and denote 

xi = vt + v(o + b)t ■ xi + xi , . . . , x ni = vt + v(o + b)t ■ x\ + x ni , 
x ni +i =vt+ yj{a + b)t- x 2 + x ni+1 , ...,x n2 =vt+ y/{a + b)t ■ x 2 + x n2 , 

Xnh-x+i =vt+ y/(a + b)t- x h + x nh _ 1+ i, ...,x nh =vt+ yj{a + b)t- x h + x nh . 
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Then 

lim ((a + b)t) h/2 P t { Xl , . . . ,x N \ yi ,...,y N ) 

= II ( 1 -<T 1 )^"* 4+l+ * i -det[«*-i(*j)]fc,J=i,...,fc J 
ie{i,...,JV}YD 

and i/ie convergence is uniform for X\, . . . ,Xh belonging to compact sets. 

Proof. Our starting point is an explicit formula for Pt(x\, . . . , a; at | Hi, ■ ■ ■ , Un) 
from [4]; it is a generalization of a similar formula for TASEP from (41]. We 
have 

P t (a;i,...,a;iv | yi,-..,2/jv) 

= ^n^"" 1 ex P (-atC, - &VCi)^ det[J' fc ,,(a;, - j/ fc )]i t 7 - 1 ) 

where 

r/c-l 



27ri -/c n^=i(i-c^) 

and the integration is over a positively oriented circle Co of a small radius cen- 
tered at 0. 

Let us study the asymptotic behavior of Fk,i{x) when t — > oo and x = x{t) = 
(a - b)t + y/(a + b)t ■ x + x. 

Lemma 7.4. For t — > oo with x = x(t) — (a — b)t + iMa~+ b)t ■ x + x, we have 

F ktl (x) = exp ((a + b)t) ((a + b^y 1 ' 2 ^ 1 ' 2 (-l) n 

riie{i,...,fc-i}\z)( 1 _ CO 

x rr ~n — (® n ^ + °0-)) 

llje{i,...,i-i}\m i 

+ J2 Ckj^xK^eMbtC- 1 +at( m ), (7.2) 

me{k...l-l}\D 

where the summation is only over indices m corresponding to distinct £ m , 
C k ,i, m (t,x)=Res z ,- 1 ( 5-1 f ~ (Cm^)"" 1 exp(6t(z - C" 1 ) + at(^~ x - Cm)) 

m \nj=i(i-o*) > 

(in particular, Ck t i m (t,x) has polynomial growth (or decay) in t as t — > oo) and 

n=|{l,...,*-l}ni3|-|{l ) ... ) /-l}nD|. 
TTie remainder o(l) above is uniformly small for x belonging to compact sets. 



Remark. Observe that in the limit regime of Lemma [7. 4| we have 
|Q a exp(6t/C m + a*Cm)l = exp [t ln(C m ) + b/C m + aCm)) 



exp t((6 - a) ln(Cm) + &/ Cm + aC, m ) 
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Since the function (6 — a) ln(r) + b/r + ar of r G R+ has a minimum at r = 1, 
for any £ m ^ 1 we have (for t 1) 

lCr" exp(&VCm + ar)\ » exp((a + b)t). 

Using the fact that Ck,i. m (t, x) have polynomial growth (or decay) in t as t — > oo 



we conclude that in the asymptotic decomposition (7.2 1 the first term is small 
comparing to the other ones. 



Proof of Lemma \7.4\ First, suppose that b > 0. Let us deform the integration 
contour so that it passes near the point 1. Since Q > 1, we need to add some 
residues: 



F k i{x) = -L f H|=i Q — ^-i exp(6tz + atz- v )dz 
2m Jci Uj^-Qz) 



+ £ 

m£{fe,...,(-l}\D 



J 

Res 



( B^-in r\ zX " eMbtz + atz ~ 1] ) ' (7 - 3) 

VIL=i(i-G*) / 



where the contour Ci is shown in Figure [4] Here the summation is only over 
indices m corresponding to distinct £ m . 





Y 







Figure 4: Contour of integration C\. The radius of the small arc is t 1 / 2 and 
the radius of larger arc is a parameter R that we choose later. 



Observe that if I < m < k and ( m > 1 , then 
Res 



= C k .i tm (t,x)C^ a; exp(6tC m 1 +at( m ). 
In particular, if the pole at C" 1 is simple (which is true, for instance, if all 
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Cm 7^ 1 are mutually distinct), then 

Res ( n,ti(i-c^) z x-i , Mg + atz ~u 



Observe that 

z*' 1 expibtz + atz' 1 ) = exp(t(\nz(a-b+r 1/2 xy / a + 6 + i -1 (x-l)) + bz+a/ ' z)) 
and 

5R((a-6) In z + 6z + a/z) 

has a saddle point at z = 1. We claim that for large t, only a small neigh- 
borhood of z = 1 gives a non-negligible contribution the integral. Indeed, let 
z = ii + iv. Then 

5R((a - 6)lnz + bz + a/z) = (a - b) In \J u 2 + v 2 + bu + (7.4) 

\/u 2 + v 2 

Setting u = 1 and differentiating with respect to v we see that for a small enough 
e > the function 

(a - b) In \J\ + v 2 + b+ -. = 

Vl + v 2 

increases on [— e, 0], decreases on [0, e] and its maximum is a + b. 

Now set the radius of the bigger arc R to be equal to + e 2 and note that 
along this arc 

(a - b) In \Ju 2 + v 2 + bu + = (a - b) In R + bu + -— 

v« 2 + v 2 R 

< (a - b) In R + b + — = (a - b) In V 1 + £ 2 + 6 + 



Summing up, if we choose any 6 > 0, then everywhere outside the 5- 
neighborhood of z — 1 the function /(z) = 3?((a — 6) lnz + bz + a/z) is smaller 
than /(l + 6) which, in turn, is smaller than /(l) = a + b. Therefore, the in- 
tegral along our contour outside the (S-neighborhood of z = 1 can be bounded 
by const ■ exp((a + b)t) exp(i(/(l) — /(l + 8))) and, thus, as t — > oo it be- 
comes exponentially smaller than exp((a + b)t), which is the smallest term in 
our asymptotic expansion (|7.2|. Consequently, this integral contributes only to 



o(l) term in (7.2 1 and can be neglected. 

To calculate the integral in the small ^-neighborhood of 1 we change the 
integration variable z = 1 + t~ x l 2 u and arrive at the integral 

t _ 1/2 j_ / naH^+£!N) 
2^4 n5iia-o(i+*- 1/a «)) 

x exp^( ln(l + t~ 1/2 u)(a -b + r l/2 Z^a + b + t~ x (x - 1)) 

+ 6(1 + r 1/2 u) + a/(l + r 1/2 u))) du 
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with contour C 2 that is a part of contour C2 shown in Figure [3] between points 
u = ±iSt 1 / 2 

Simplifying the integrand we arrive at 

rV«p W . + t )) ; 1 I nt.'d-Od + r" 2 ")) 



2ttz 



c 2 n5=i(i-<7-(i+*- l/2 «)) 



x exp(u 2 (a + b)/2 + uxVa + b + o(l)) . 



Making the linear change of variables v = u\/a + b and sending t to 00 we obtain 
the required integral. 

As for the case 6 = 0, the argument is similar and we omit it. The only 
difference is that now contour of integration should have no vertical line part, 
i.e. it consists of the arcs of unit circle and i _1 / 2 -circle. □ 



Now we continue the proof of Proposition 7.3 

The next step is to do certain elementary transformations to the matrix 
\F^i{xi — yk)] in order to simplify its determinant. First, suppose that all Cm 
(except for those equal to 1) are distinct. 

We take the (N — l)st row of matrix [Fk.i(xi — yk))ki=i (i- e - k = N — 1), 
multiply it by 

rv n -yN-i \ (7 K\ 

llj=n I 1 _ Ci/Gv-iJ 



and add to the nth row for n = 1, . . . , N—2. As a result the term in (7.2 1 coming 
from the residue at Cjv-i remains only in the (N — 1, N) matrix element. 

Next, we take the (N — 2)nd row of the matrix and add this row (again 
with coefficients) to rows 1, . . . , N — 3 so that the term in ( |7.2| coming from 
the residue at remains only in the (JV — 2, JV — 1) and (N — 2, iV) matrix 

elements. 

Repeating this procedure for every row h such that Ch ^ 1 we get a trans- 
formed matrix [Fkj(xi ~yk)]ki—i- I n this new matrix, the term in (7.2 1 coming 



from the residue at £/i remains only in the matrix elements Fh-i t h{xh ~Vh~ 1); 

■ ■ ■ , Fh-\,n{xn - Vh-i)- 

Let us now do some columns transformations. We take the second column 
(1 — 2) and add it with coefficients to the columns 3, . . . , N so that the terms 
coming from the residue at £i vanish in columns 3, . . . , N. Then we repeat this 
procedure for the third column and so on. 

Finally, we get amatrix [Fk t i(xi — yk)]k i=\ whose determinant coincides with 
that of [Fk,i(xi — yk)]k i— it but this matrix has only N — h elements with order 
of growth greater than exp((a + b)t) (i.e. elements with terms coming from the 
residues). These elements are F n ^i_ n (x n — y n -i) for n € {2, . . . , N}, n £ D + l, 

F n - hn (x n - y n _i) = -Q^Cn'-T^- 1 exp(6tC»_i + ^C«-i)(l + o(l)). 

It follows that asymptotically as t —¥ 00 the determinant det[Fk t i(xi — yk)\ 
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factorizes: 



det[F k ,i(xi - yk)]k,l=i,...,N 

IJ C-l +Vn " eMbK- 1 + atCm) det[G fc ,,] fc)J= i,... lh (l + o(l)), 



ne{2,...,JV}\(D+l) 



(7.6) 



where [Gfc^Jfc^i,...,/, is the submatrix of [Fk,i(xi - yk)] with rows ni,n 2 , ■■ .,n h 
and columns 1, rii + 2, ri2 + 1, . . . , Uh—i + 1. Looking back at row and column 
operations that we made with [Fkj(xi — yk)] to get [Fk,i{xi — yk)], we see that 
the determinant of Gk,i coincides with the determinant of the Gaussian terms of 



the decomposition (7.2) of matrix elements Fk t i(xi — yk) in rows ni,n%, . . . ,rih 



and columns 1, n\ + 1, n 2 + 1, . . . , Uh—i + 1. Therefore, 

det[G fc>I ] fcl , = i,..., fc = exp(M(a + b))(-l) h ((a + b)t)- h ' 2 (o(l) 

h h-1 

+ u n a- con n det^^o]^!,.^ 

m=lie{l,...,n m -l}\D m=l jG{l,...,n ro }\-D ' 

= cxp(ht(a + b))(-l) h ((a + b)t)- h / 2 l[ (1-Ci) 

ie{i,...,JV-i}\£) 

x (det[0 fc _i(xi)]fc,l=i,-,fc + (!)) • ( 7 - 7 ) 



Combining ( 7.1 1 with asymptotic formulas ( |7.6| and ( |7.7| we conclude the proof. 



Now suppose that some of Q coincide. In this case the computation of the 
residues in the decomposition ( |7.2[ ) becomes more complicated, since some of the 
poles are not simple. However, the scheme of the proof remains the same: First 
we transform the matrix [Fk t i(xi—yk)] by means of elementary row transforms so 



that the terms in (7.2) coming from the residue at remain only in the matrix 
elements F h _ lth (x h - y h - 1), . . . , F h _i tN (x N - y h -i)- Then we do elementary 
column transforms so that only N — h elements with order of growth greater 
than exp((a + b)t) remain in the resulting matrix. After that the computation 
of the determinant repeats the case of distinct Q m . 

□ 
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